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ABSTRACT 


The  point  availability  of  a  system  of  components  each  of 
which  is  subject  to  random  failures  and  has  random  restoration 
times  is  determined.  Each  component  is  assumed  to  have  a 
fixed  number  of  spares  such  that  when  all  spares  are  exhausted 
no  restoration  can  take  place. 

Exact  expressions  are  obtained  for  the  availability  of  any 
series/parallel  system  composed  of  independently  operating 
components.  Furthermore,  expressions  are  obtained  for  the 
availability  of  series  systems  under  five  operating  scenarios 
for  which  components  sure  not  independent.  Order  relationships 
on  system  availability  for  the  various  operating  scenarios  are 
determined. 

An  algorithm  which  optimally  allocates  spare  parts  over  the 
components  in  a  system  is  presented.  The  algorithm  offers 
computational  savings  over  solutions  derived  from  dynamic 
programming . 
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LIST  OF  SYMBOLS 


A: 
n : 

F  ( t)  : 

f  (t)  : 
G(t)  : 

g  (t) 

£(n)c>  = 

f  *g: 

A: 


D 


2 


B: 


r*  (s) 


Constant  failure  rate  of  component  parts 
Constant  replacement  rate  of  spare  parts 
Distribution  function  of  lifetime  of  component 

parts  (l-e"^t) 

dF(t)  _  -At 
^t  Ae 

Distribution  function  of  replacement  time  of  spare 

parts  (l-e_nt) 

dG ( t)  _  _-nt 

"dt—  -  ne 

n-fold  convolution  of  f 
Convolution  of  f  and  g 

Two  components  are  in  "up"  states  in  a  two-component 
series  system  when  each  component  forms  an  alternating 
renewal  process  ("up”  and  "down") 

Component  #1  is  in  "down"  state  and  component  #2 
is  in  "up"  state  in  the  above  case  (see  "A") 

Component  #2  is  in  "down"  state  and  component  #1 
is  in  "up"  state  in  the  above  case  (see  "A") 

Both  components  are  in  "down"  states  in  the  above 
case  (see  "A") 

/  e"tsr(t)dt 
0 


Limiting  probability  of  state  i  in  a  Markov  Chain 
P^:  Transition  probability  matrix  of  an  embedded  Markov 

Chain  under  scenario  i 
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q^j/k/S)  = 


[qi(2n)] 


p»q 


A(n) 


(t)  : 


prob{X„  =  j /X  .  =  k,t  -  t  =  s)  for  scenario 
n  n+l  n+i  n 

i  where  xm  is  state  of  embedded  Markov  Chain  after 
m 

mth  transition  and  t  is  point  in  time  of  mth 

m 

transition 

Element  (p,q)  of  2n-fold  convolution  of  matrix 
*i 

Availability  of  one-component  system  when  there 


are  n  spares  available  starting  from  state  "up" 
at  time  0 


A<n)  (t) 
av 


B(n|(t) 


A<nl 


(t)  : 


a(n,m) 

Ai 


(t)  : 


A|i'k>Dj (t) : 


Average  availability  A^  (t)  =  4  /  A^  (s)ds 

av  z  o 

Availability  of  one-component  system  when  there 
are  n  spares  available  after  time  0  starting  from 
"down"  state  with  original  item  failed  on  or 
before  0 

Availability  of  a  system  at  time  t  under  scenario 
i  with  every  component  started  from  "up"  state 
at  time  0  when  there  are  n  spares  usable  for 
every  component  of  the  system 
Availability  of  a  two-component  series  system 
at  time  t  under  scenario  i  with  both  components 
started  from  “up"  states  at  time  0  when  there 
are  n  spares  available  for  component  #1  and  m 
spares  available  for  component  #2 
Availability  of  a  two-component  series  system  at 
time  t  under  scenario  i  with  the  process  started 
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from  state  D  ^  or  time  0  when  l  spares  are 
available  for  component  #1  and  k  spares  are 
available  for  component  #2 

A±(t)  =  lim  a!"1 (t>  =  lim  A|n,m) (t) 

n-+°°  n-*-°° 

m-+°° 

A.  =  lim  A. (t) 

A ( 2 ,m' n) ( t) :  Availability  of  three  component  system  (for  any 
configuration)  at  time  t  with  every  component 
of  the  system  started  from  "up"  state  at  time 
0  when  there  are  i  spares  for  component  #1,  m 
spares  for  component  #2  and  n  spares  for 
component  #3 

A(n,m) .  Availability  of  two  component  series  system  at 
s 

time  t  with  both  components  started  from  "up" 
states  when  two  components  are  mutually  indepen¬ 
dent  and  there  are  n  spares  available  for 
component  #1  and  m  spares  are  available  for 
component  #2 

•A^n,m) (t) :  Availability  of  two  component  parallel  system 

at  time  t  in  the  above  case  (see  (t) ) 

s 

a|n^ (t) :  Stage  i  availability  (stage  return)  in  dynamic 

programming 

Ai^i^:  Maximum  i-stage  return  in  dynamic  programming 

s  t 

Xi  ^  X2:  Random  variable  is  stochastically  larger 

than  random  variable  X2 
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INTRODUCTION 


I . 

A  key  concern  of  military  planners  is  the  material  readi¬ 
ness  of  their  weapon  systems .  Many  attempts  have  been  made 
in  the  past  to  develop  measurable  indicators  of  material 
readiness.  Among  these  are  reliability,  fill  rate,  down  time, 
maintenance  manhours,  supply  response  time,  and  availability. 

Out  of  these  attempts  to  measure  readiness,  operational  availa¬ 
bility  has  become  the  most  widely  accepted  indicator.  Ref. 

[32]  states,  "while  it  is  an  indicator  that  operating  forces 
can  recognize  and  relate  to,  it  also  has  the  feature  of  reflecting 
in  a  logical  manner  the  relationship  of  the  different  factors 
which  contribute  to  material  readiness."  Much  effort  is  under¬ 
way  trying  to  determine  how  weapon  systems  should  be  built, 
maintained,  and  supported  with  spares  in  order  to  maximize 
availability  with  given  resources.  This  thesis  examines  some 
of  the  problems  concerned  with  determining  the  operational 
availability  of  a  weapon  system  and  then  allocating  scarce 
resources  to  the  system  to  optimize  availability. 

Many  people  talk  and  write  about  operational  availability, 
but  not  all  have  the  same  notion  as  to  what  is  meant  by  the 
expression.  Operational  availability  is  generally  considered 
to  be  a  measurd  of  the  likelihood  that  a  system,  when  used 
under  stated  conditions  in  the  actual  operational  environment, 
shall  operate  satisfactorily.  It  is  clear  that  component 
reliabilities,  maintenance  policy,  spare  parts  support,  system 
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configuration,  repair  times,  and  supply  response  times  all 
have  impact  on  availability.  However,  this  notion  of  operational 
availability  can  give  rise  to  more  than  one  workable  mathemati¬ 
cal  expression.  To  illustrate,  consider  a  ship  on  a  60-day 
deployment.  Let  A(t)  be  the  probability  that  a  given  weapon 
system  is  operational  at  time  t.  One  might  use  any  of  the 
following  as  reasonable  measures  of  availability: 

a)  A (60)  (end  of  deployment  point  availability) 

x  60 

b)  Jo  /  A(t) dt  (average  availability) 

c)  The  probability  that  the  system  shall  operate 
satisfactorily  when  called  upon 

,,  Mean  time  to  failure  (MTBF) _ 

'  MTBF  +  Mean  time  to  replace  (MTTR) 

The  first  expression  is  the  likelihood  that  the  equipment  will 
be  operational  at  the  end  of  the  deployment;  the  second  repre¬ 
sents  an  average  availability  over  the  duration  of  the  deploy¬ 
ment;  the  third  considers  oniy  the  likelihood  that  the  system 
is  operational  at  those  times  when  needed  (thus  introducing 
the  mission  duty  cycle  as  a  factor  to  consider) ;  and  the  last 
expression  is  the  "definition"  of  availability  specified  by 
Ref.  t)2] ,  and  used  most  frequently  in  practice.  Under  certain 
very  restrictive  conditions  the  above  expressions  may  yield 
essentially  the  same  values.  However,  for  most  real-world 
cases  the  expressions  are  not  equivalent;  they  can  be  sub¬ 
stantially  different.  The  following  specific  cases  illustrate 
some  differences. 
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Case  1:  The  equipment  is  a  single  component  which  has 


exponential  life  time  with  failure  rate  x .  There  are  unlimited 
spares  aboard  ship.  The  distribution  of  the  sum  of  supply 
response  time  (time  required  to  get  the  replacement  part  to 
the  equipment)  and  repair  time  (hereafter  called  replacement 
time)  is  exponential  with  rate  n.  The  system  is  operational 
at  the  start  of  the  deployment. 

Define  the  random  process  (X(t)  :  t  ^  0}  as  follows: 


X(t) 


'  1  if  component  is  operational  at  t 

< 

.  0  otherwise 


Under  the  conditions  above,  X(t)  is  a  renewal  process  and 
A(t)  =  P{X  ( t)  =1]  is  easily  derived  from  renewal  theory  to 
be 


A(t)  =  -r-2—  +  exp{-  (X+n)  t} 

A  +  n  A  +  ri 


(1-1) 


or,  since  MTBF  =  i  and  MSRT  +  MTTR  =  —  we  can  write  (1-1)  as 

x  n 


A  (t ) 


MTBF 

MTBF  +  (MSRT  +  MTTR) 


+ 


(MSRT  +  MTTR) 
MTBF  +  (MSRT  +"  MTTR) 


exp{- (x+n) t} 


(1-2) 


To  allow  numerical  comparisons  of  the  four  interpretations  a) , 
b) ,  c) ,  and  d)  of  availability,  let  x  =  0.01  and  n  =  0.10.  Then 
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a)  A{60) 


0.90921 


60 


■»  u  u 

b)  ^  /  A(t) dt  =  0.92285 


MTBF 


d)  MTBF  +  (MSRT  +  MTTR)  =  °*90909 


There  are  only  minor  differences  between  a)  and  d)  but  a  bit 
larger  difference  between  b)  and  the  others.  The  calculation 
for  interpretation  c)  yields  exactly  the  same  result  as  b)  if 
the  times  that  the  component  is  needed  are  uniform  on  the  inter¬ 
val  (0,60]  .  Interpretation  d)  provides  the  most  conservative 
estimate  of  availability.  From  equation  1-2  one  can  easily 
see  that 


lim  A ( t) 


MTBF 

MTBF  +  (MSRT  +  MTTR) 


Thus  d)  is  simply  the  limiting  availability .  This,  however, 
is  not  always  the  case  as  is  shown  later. 

A  plot  of  A(t)  vs.  t  demonstrates  how  point  availability 
vaies  as  a  function  of  time  (length  of  deployment) .  See 
figure  1.1. 
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Consider  now  a  case  that  is  more  operationally  realistic. 
Significant  differences  will  be  seen  in  the  estimates  yielded 
by  the  different  expressions. 

Case  2:  The  equipment  is  a  single  component  which  has 
exponential  life  time  with  rate  x .  There  is  a  single  spare 
aboard  ship  and  no  chance  for  resupply  until  the  deployment 
is  over.  The  replacement  time  is  exponential  with  rate  n. 

The  system  is  operational  at  time  0. 

The  random  process  {X(t),t  >  0}  which  describes  the  up/ 
down  status  of  the  equipment  is  not  now  a  renewal  process. 

A  sample  path  for  X(t)  is  shown  in  figure  1.2. 


Fig.  1.2  Sample  path  for  X(t)  when  only  one  spare  is 
stocked 


The  important  thing  to  note  is  that  the  component  will  remain 
in  a  down  status  after  time  t^  since  there  is  no  spare  for 
replacement.  In  the  next  chapter  an  expression  is  derived 
for  the  point  availability  from  which  the  following  values 
are  determined  (x  =  0.01  and  n  =  0.10): 
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a)  A(60) 


0.8472 


60 


b)  ^  /  A ( t) dt 


d) 


0 

MTBF 


MTBF  +  MTTR 


=  0.9054 

0.9091 


Greater  differences  among  the  availability  estimates  are 
observed  when  a  single  spare  is  stocked.  The  differences 
can  be  made  quite  significant  by  appropriate  choice  of  the 
parameters  X  and  n.  For  example,  if  X  =  0.05  and  n  =  0.50, 

a)  A{60)  =  0.2096 

.  60 

b)  /  A(t)  dt  =  0.578361 

j \  MTBF  _  anQ1 

d)  MTBF  +  MTTR  ~  0.9091 


The  above  examples  show  that  there  can  be  significant 
differences  in  the  values  obtained  from  the  different  defini¬ 
tions  of  availability.  Each  expression  has  its  supporters. 
For  logistics  planning  purposes,  military  planners  want  to 
determine  the  likelihood  that  a  deployed  unit  will  complete 
a  mission  with  a  given  system  operational.  Therefore,  they 
need  to  know  the  point  availability  at  time  t  where  x  is  the 
length  of  the  deployment  period.  Furthermore,  the  average 
availability  and  the  limiting  availability  acre  both  functions 
of  the  point  availability.  Therefore,  we  focus  primarily  on 
point  availability  in  this  thesis.  Examples  of  calculations 
of  average  availability  and  the  limiting  availability  are 
provided  to  illustrate  their  calculation  .  The  definition  of 
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availability  based  on  the  ratio  of  mean  time  to  failure  to 
the  sum  of  mean  time  to  failure  and  mean  time  to  replace 
leaves  much  to  be  desired.  The  definition  assumes  that  the 
up  and  down  states  of  the  system  satisfy  an  alternating 
renewal  process.  This  implicitly  assumes  infinitely  many 
spares  making  the  expression  useless  for  determining  the  number 
of  spares  required  to  support  a  system.  Furthermore,  the  ratio 
expression  is  a  limiting  result,  whereas,  the  major  interest 
is  in  the  availability  at  specified  points  in  time  or  over 
specified  intervals.  Finally,  the  expression  does  not  lend 
itself  for  calculating  the  system  availability  as  a  function 
of  the  availabilities  of  its  components. 

Let  us  now  look  at  the  point  availabilities  of  some  simple 
systems.  Consider  first  the  case  in  which  the  system  is 
composed  of  two  components  connected  in  series  each  with 
infinitely  many  spares.  We  assume  that  the  components  operate 
independently.  It  is  easy  to  see  from  first  principles  that 
the  system  availability  is  the  product  of  the  component 
availabilities.  That  is, 

Asys(tl  "  Al(t)  '  A2(t) 

Similarly,  if  the  system  is  composed  of  two  components  in 
parallel,  the  system  availability  can  easily  be  shown  to  be 

Asys(t)  =  1  -  U-AjCt))  (1-A2(t))  . 
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The  reader  familiar  with  reliability  theory  will  recognize 
the  above  calculations  as  identical  to  those  used  to  calculate 


reliabil’ti  of  systems  composed  of  two  components  in  series 
and  parallel,  respectively.  These  results  are  easily  extended 
to  any  finite  number  of  components.  The  formulae  are  also 
valid  when  each  component  has  only  finitely  many  spares  for 
support.  Thus,  the  system  availability  of  any  series/parallel/ 
mixed  system  can  be  determined  easily  from  the  availabilities 
of  its  components  provided  the  components  operate  independently. 

Whenever  the  components  of  a  system  do  not  operate  inde¬ 
pendently  the  formulae  given  above  are  not  correct.  There 
are  many  real  world  examples  in  which  the  operation  of  one 
component  may  depend  on  the  operation  of  another  component. 

When  this  happens,  the  calculations  are  much  more  complicated. 
First,  the  nature  of  the  dependence  must  be  described.  Then, 
conditioned  on  the  nature  of  the  dependence,  the  component 
availabilities  must  be  aggregated  appropriately  into  an 
expression  for  system  availability.  Certainly,  one  cannot 
analyze  all  possible  types  of  dependence.  However,  there  are 
certain  types  which  are  most  realistic  from  operational  considera¬ 
tions.  In  this  thesis,  we  consider  explicitly  four  cases, 
called  operational  scenarios,  in  which  the  component  operations 
are  interrelated.  The  examples  below  demonstrate  the  differ¬ 
ences  in  system  availability  caused  by  a  change  in  scenario. 

The  formulae  from  which  the  calculations  are  made  are  derived 
in  this  thesis. 
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Case  3: 


Consider  a  two-component  series  system  with  both  com¬ 
ponents  identical.  We  assume  that  the  component  lifetimes 
and  replacement  times  are  independent  and  exponentially 
distributed  with  parameters  X  =  1/30  and  n  =  1/5,  respectively. 

We  consider  first  the  case  in  which  the  components  operate 

<n.,n2) 

independently.  Let  A  (t)  be  the  system  availability  at 

time  t  when  there  are  n^  spares  for  component  1  and  n2  spares 
for  component  2.  We  determine  in  this  thesis  the  following 
availabilities  for  t  =  90. 


Case  4: 

Consider  now  the  same  system  configuration,  component 
reliabilities,  and  component  replacement  times.  However, 
now  suppose  that  the  surviving  component  shuts  down  whenever 
a  failure  occurs.  Under  this  operational  scenario,  the  following 
availabilities  are  obtained: 
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(n1 ,n0) 


<nl'n2} 

A  1  ^  (90) 

(0,0) 

0.0025 

(1,0) 

0.013 

(1,1) 

0.059 

(2,1) 

0.129 

(2,2) 

0.266 

We  observe  that  the  availabilities  are  slig..  -eater  than 
in  the  previous  case  where  the  components  operate  independently. 

Case  5; 

For  the  last  case,  consider  the  same  situation  described 
above,  but  assume  that  the  identical  components  share  a  common 
pool  of  spare  parts.  Let  A^(t)  be  the  system  availability 
at  time  t  when  there  are  n  spares  in  the  pool.  The  following 
values  are  determined  in  this  thesis: 


n 

A  n  (t) 

0 

0.0025 

1 

0.023 

2 

0.094 

3 

0.234 

4 

0.412 

Given  a  fixed  number  of  spares,  we  observe  that  system  availa¬ 
bility  is  greater  when  spares  are  shared  than  when  each  com¬ 
ponent  has  unique  support. 
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The  examples  and  discussion  above  illustrate  that  the 
availability  of  a  system  composed  of  multiple  components 
arranged  in  series/parallel/mixed  configurations  depends  on 
many  factors.  Among  these  are: 

(1)  component  reliabilities 

(2)  component  replacement  times 

(3)  system  configuration 

(4)  operational  scenario 

(5)  number  of  spares 

(6)  spares  sharing  policy 

Any  effort  to  calculate  system  availability  that  does  not 
explicitly  consider  each  of  these  factors  cannot  be  accurate. 

The  expressions  derived  in  this  thesis  take  into  account  each 
of  these  factors.  As  a  result,  the  expressions  derived  are 
complicated,  but  that  is  frequently  the  price  paid  for  accuracy. 
Only  after  the  exact  expressions  are  known  can  any  serious 
effort  be  made  to  obtain  reasonable  approximations  and  simpli¬ 
fications  . 

In  the  remainder  of  this  chapter  we  describe  the  results 
obtained  in  the  thesis.  For  mathematical  simplicity  and  con¬ 
formance  to  various  military  specifications  we  assume  that  both 
the  lifetimes  and  the  replacement  times  of  the  individual  com¬ 
ponents  are  independent  and  identically  distributed  (iid)  with 
exponential  distributions. 
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In  Chapter  II,  computational  expressions  for  point  avail¬ 
ability  and  average  availability  of  a  one-component  system 
having  a  fixed  number  of  spares  are  determined. 

The  tradeoff  of  repair  versus  supply  is  examined  by  con¬ 
sidering  two  policies;  1)  a  replacement  policy  with  a  fixed 
number,  n,  of  spares  and  2)  a  repair  policy  with  repair  rate 
n'  (implicitly  equivalent  to  an  assumption  of  infinitely  many 
spares) .  One  aspect  of  this  tradeoff  is  discussed  by  looking 
at  what  the  repair  rate  n ’  must  be  in  order  to  achieve  the 
same  level  of  availability  provided  by  the  replacement  policy 
with  n  spares . 

Given  a  fixed  length  of  mission  duration  and  a  finite 
number  of  spares,  a  system  may  not  be  available  at  the  end  of 
a  mission  due  to  lack  of  spares.  The  probability  distribution 
of  system  down  time  due  to  lack  of  spares  is  determined  as  a 
function  of  number  of  spares  and  mission  duration. 

In  Chapter  III,  we  examine  two-component  series  systems 
under  various  operational  deployment  scenarios.  We  consider 
the  case  in  which  the  two  components  operate  independently 
nad  four  cases  in  which  the  components  do  not  operate  inde¬ 
pendently.  The  four  dependent  scenarios  are  distinguished 
by  what  happens  to  the  surviving  component  when  a  failure 
occurs.  The  scenarios  considered  are  all  operationally 
realistic.  We  consider  cases  in  which  there  are  infinitely 
many  spares  and  finitely  many  spares  for  each  component.  We 
also  consider  the  case  in  which  there  are  finitely  many  spares 
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which  are  shared  by  the  two  components.  Exact  results  are 
obtained  for  the  system  point  availability  in  some  cases, 
and  approximate  results  in  others. 

In  Chapter  IV,  exact  computational  formulae  are  derived 
for  the  availability  of  series/parallel/mixed  system  composed 
of  independently  operating  components  when  each  component  has 
a  fixed  number  of  spares.  In  addition  an  approximation  is 
obtained  for  the  availability  of  a  two-component  series  system 
operating  under  the  "symmetric  shut  down"  scenario  when  each 
component  has  a  unique  finite  spares  pool. 

In  Chapter  V,  order  relationships  are  established  on  sys¬ 
tem  availability  for  the  five  operating  scenarios.  The  exact 
computational  expressions  for  the  availability  of  a  system 
under  the  non-independent  scenarios  are  very  complex.  How¬ 
ever,  these  order  relationships  can  be  used  for  establishing 
approximations  or  bounds  for  the  exact  system  availability. 

In  Chapter  VI,  the  results  derived  previously  are  utilized 
in  developing  an  algorithm  for  optimally  allocating  a  given 
budget  for  spare  parts  support  for  the  components  of  a  system. 
The  algorithm  views  the  system  as  a  series/parallel  mixture 
of  independently  operating  modules.  A  module  can  consist  of 
components  which  do  not  operate  independently.  The  algorithm 
will  handle  any  system  for  which  the  dependently  operating 
components  are  described  by  the  scenarios  addressed  in  this 
thesis.  The  algorithm  is  developed  from  dynamic  programming. 
The  result,  however,  is  an  algorithm  that  offers  substantial 
computational  efficiency  over  dynamic  programming. 
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The  result  is  a  procedure  that  could  be  utilized  in  the 
construction  of  allowance  or  load  lists.  A  comprehensive 
example  is  presented  and  solved  using  the  algorithm. 

Chapter  VII  summarizes  the  results  of  this  research 
effort  and  discusses  problems  of  implementation.  It  recommends 
additional  work  for  subsequent  research  efforts. 
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II.  AVAILABILITY  OF  ONE  COMPONENT  SYSTEM  WITH 
FINITE  NUMBER  OF  SPARES 

A.  POINT  AVAILABILITY 

In  this  section  we  determine  mathematical  formulae  for  the 
point  availability  of  a  single  component  system  having  n  spares 
and  no  repair  capability.  Since  there  are  only  finitely  many 
spares  and  no  possibility  for  repair,  the  system  will  alternate 
between  up  and  down  states  until  all  spares  are  exhausted  and 
will  remain  in  a  down  state  whenever  the  last  spare  fails. 

Let  us  introduce  some  notation  that  will  be  used.  Let  T^ 
be  the  lifetime  and  the  replacement  time  of  the  ith  unit. 

We  assume  that  {T^}?_^  are  independent  and  identically  distributed 
(iid)  with  distribution  function  F(t)  and  probability  density 
f(t).  Similarly,  the  replacement  times  are  iid  with 

distribution  G(t)  and  density  g(t).  Furthermore,  the  replace¬ 
ment  times  are  independent  of  the  lifetimes.  We  use  f*g  to 

( k) 

indicate  the  convolution  of  f  and  g  and  f  to  represent  the 
k-fold  convolution.  Finally,  let  Pk(t)  be  the  probability  that 
the  kth  unit  will  be  in  operation  at  time  t. 

The  first  result  we  derive  is  the  general  expression  for 
the  point  availability  of  the  single-component  system  when 
there  are  n  spares,  A^  (t)  . 

The  system  will  be  operational  at  time  t  if  and  only  if 
the  kth  unit  (k  =  1,  2,...,n+l)  is  in  operation  at  time  t.  Thus, 

/ _ \  n+1 

A  n  (t)  =  l  pk(t)  . 

k=l 
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Now 


k-1  k-1 

P.(t)  =  P[  l  (T.  +  R • )  <  t  and  T.  +  l  (T.  +  R. )  >  t] 

k  i=l  1  1  K  i=l  1  1 

=  f  (i71)  ?(t_s)  ds 

=  [ (f *g) (k'1)*P] (fc) 

where 

F(t)  =  >  t]  . 

We  have, 

A(n)(t)  =  l  [ (f*g) (k_1)*F)  . 

k=l 

or 

A (n)  (t )  =  F(t)  +  l  [ (f*g)  (k)*F]  .  .  (2-1) 

k=l  m 

This  result  can  be  rewritten  recursively  as 

A(n) (t)  =  A(n_1) (t)  +  [(f*g) <n)*F](t)  (2-2) 

and 

A (n) (t )  =  F (t)  +  [A(n_1)*f *g] (fc)  (2-3) 
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Each  expression  has  some  usefulness.  Equations  (2-1)  and  (2-2) 
are  probably  preferable  computationally.  Equation  (2-2)  pro¬ 
vides  a  simple  expression  for  the  marginal  contribution  that 
the  nth  spare  provides  to  system  availability. 


A (n) (t)  =  A(n) (t)  -  A(n_1) (t)  =  t(f*g)(n)*F] 


(t) 


For  the  special  case  in  which  f(t)  =  Xe  and  g(t)  =  ne~nt 
the  convolution  expression  is  found  using  Laplace  transforms 
to  be 


(f  ^  *g  ^  *f] 


,8.k  rtk  k  .  ,.r  (k+r-l)Pr  t*k  ,  -Xt 
(t)  -  +  TE^iTT1* 


(k-1)  k-1  (k+l)pt  ,  ■ -nt 


1,  a  ^  ^  Vrw— j-/  _  J-  -mu  /0  AS 

<2'4) 


where 


nPk  =  ■,  ,  6  =  x  n ,  and  5  =  n  -  X  >  o . 

As  an  example,  if  x  =  1/30,  n  =  1/5  and  t  =  90,  we  have 


[  ( f *g) (7,*F] 


(t) 


rei 7 ,t7  7  t6  7-8. t5  7 -8  •  9  .  t4 

lV  [TT  7  FT  777?  TT  "YTfiT  77 


2 1 6  * 


7  ’  8  •  9  *  10  t3  7  •  8 • 9  •  10 • 11  t2  ,  7  •  8  •  9  •  10 • 11 *12  t 


416 


3  ! 


516 


616 


*17 


7 -8 • • -12  *13, _-X t  .  ,e7, .  ,t6  ,  8  t5  .  8-9  t4  .  8-9-10  t3 

— 7^7 - +  [t§-]  lFr  +  T,Tr  +  77^‘47  +  _3!6. 


T ~'TT 


8-9-10-11  8-9---11 

.4  *21  + 


416 


- ^Ai.t  *  i2--3]«~nt  -  0.0033 

5 1 6  3  616° 
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as  the  marginal  contribution  of  the  7th  spare  to  the  availa¬ 
bility  at  t  =  90. 

An  observation  of  the  convolution  expression  provides  us 
with  a  useful  approximation . 

[{f*g) (n)*F] (t)  =  I[fn+1*gn] (t) . 


Now  the  (n+l)-fold  convolution  of  f  is  the  gamma  density  with 
parameters  (n+l,i),  and  the  n-fold  convolution  of  g  is  the 
gamma  density  with  parameters  (n,  i  )  .  As  n  gets  large  these 

converge  to  normal  densities  with  parameters  and 

n  n  A  x  1 

(— ,  respectively.  Thus,  for  large  n,  the  convolution 

n  1 
expression  be  be  approximated  by  — h(t),  where  h  is  a  normal 

density  with  mean  +  —  and  variance  +  JL  .  in  terms 

x  n 

of  the  standard  normal  density  $ ( • )  we  have 


[(f*g) (n)*F] (t) 


/n2(n+l)+\2n 


t-  (2fi  +  a) 

♦  [  ■  -— * - S.-] 


\  n 

For  the  example  above,  with  n  =  7,  \  =  1/30,  =  1/5  and  t 

we  have 


=  90, 


A(7){90)  =  0.013 

Compared  to  the  exact  value  of  0.0033  the  approximation  is 
not  very  good,  but  it  would  improve  as  n  gets  larger. 
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B.  AVERAGE  AVAILABILITY 

Let  A^(t)  be  the  average  availability,  i.e., 

av 


*av  «=> 


-  I  A(n) (s)ds 

*•  A 


Then,  from  eq.  (2-1) ,  (2-2) ,  (2-3)  and  the  identity 


/  t“e  d 


m  _ .  m-p 

_ax  r  /  nP  m!  t 
*  ‘  6 
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we  have 


*av  (t) 


1  n  -*t,  .  1  r  ftr0ikrxk 

It[1"e  1  +  t  VJ  [JTT 


y  .  .  r  (k+r-1) Pr  t(k‘r).  -Xt 

rh  - TT^71^ 


1  r  rfc.  ,,fc+l)  r  ek  ,  ,TW 

t  kix  y_1)  t7TT]  1  (k-TTT 


{krL)  (k+4)P4  -nt. 


+  I  k-r-TTT]e  dT 

a-1  ily  [K  '  ' 


1  n  ^-xti  j.  1  Y  reikf  r  1 

-  xt[1_e  1  +  t  {t7+T 


k  ,  ,%P  e'Xt  t(k"p>  ,  £  f  \ r .  (k+r-1)  Pr 

+  Jo1'11  +  r=l  11 


k-r  -x  t  k-r-i 

■  i  i  +  ‘kIll(-i)p-TE?^ 

fc  k=l  <5k  1  ^  p=0  lK  1 


+  <*?'  (k+*)P*f  1 

(-n)p+1  4-1  U6a  nk_iL 


k-4-1  -nt  (k-i-l-p ) 

+  I  (-DP  - - - -  T+I - ll- 

p-0  (k-i-l-p)!  (-n) p 


(2-5) 
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For  a  numerical  example,  let  A  =  yy,  n  =  y  and  t  =  90. 

From  eqs.  (2-1),  (2-2),  (2-3),  (2-4),  we  determine  the  following 
values : 


Table  2.1:  Point  Availability  vs.  Average  Availability 


\ 

A(n) (t) 

A(n) (t) 
av 

II 

o 

0.0498 

0.316737644 

1 

0.2170845 

0.57439463 

2 

0.4622845 

0.739965467 

3 

0.6721845 

0.81423 

4 

0.791011967 

0.844388 

5 

0.8388 

. 

0.8534416 

6 

0.8531 

0.8555475 

7 

0.856412 

0.856957 

00 

0.859 

0.8639455 

C.  TRADEOFF  OF  REPLACEMENT  VS.  REPAIR 

A  high  level  of  operational  availability  can  be  achieved 
in  different  ways.  One  can  increase  the  reliability  of  the 
system  components,  one  can  provide  generous  spares  support, 
one  can  build  in  system  redundancy ,  an&.one  can  provide  a 
maintenance  capability.  After  an  equipment  is  put  into  opera¬ 
tion  there  is  little  chance  to  do  anything  about  component 
reliability  or  system  configuration.  However,  one  can  still 
consider  logistics  tradeoffs  between  providing  spares  support 
and  a  repair  capability  which  consumes  only  piece-parts  support. 
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Suppose  that  a  system  can  be  repaired  with  repair  rate 
n'  and  assume  that  infinitely  many  repairs  can  be  made.  The 
repair  rate  n '  required  to  provide  a  specified  level  of  avail¬ 
ability,  A(t) ,  can  be  determined  by  solving  for  n'  in: 

V(t)  =  i+tt  +  exp_(x+n,)t  (2*5) 

In  evaluating  the  tradeoff  between  repair  and  spares  support 
it  is  useful  to  compare  the  repair  rate  n '  with  the  number 
of  spares,  n,  required  to  achieve  the  same  point  availability. 

That  is,  we  compare  o'  with  the  value  of  n  which  satisfies 

An(t)  =  exp(-At)  +  l  [(f*g) (k)*F]  ,  (2-6) 

k=l  m 

Certainly,  there  are  other  factors  to  consider  in  making  a 
decision  about  repair  vs.  replacement.  For  example,  one  would 
need  to  consider  manpower  requirements,  training,  piece-parts 
support  needs,  and  space  considerations.  However,  the  tradeoff 
discussed  above  is  useful  as  an  indicator  of  whether  or  not 
one  should  even  consider  repair. 

As  a  numerical  example  let  A  =  ,  n  =  ^-,  and  t  *  90. 

A(1) (t)  =  e"Xt  +  (f*g*F) (t)  =  0.2170845 

To  achieve  this  level  of  availability  we  find  that  the  maintenance 
rate  n*  should  be  at  least  0.008233. 
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For  n  =  2 


A(2) (t)  =  e"At  +  l  [(f*g)(k)*F]  }  =  0.4622845. 

k*l 

On  solving  (2-5)  for  n ' ,  with  Aq , (t)  set  to  0.4622845,  we 
find  that  n 1  =  0.028423. 

Tables  2.2  shows  the  comparisons  between  n  and  n'  for 
n  s  l  to  7  (for  values  of  n  larger  than  7  there  is  very  little 

increase  in  A^  (t) )  . 

Table  2.2:  Repair  Rate  as  a  Function  of  n  (for  a  given  (X,n>) 


A(t) 

Availability  U,n) 

- - 

n* 

n  =  1 

0.2170845 

0.008233  (n  =  24. 3n’) 

2 

0.4622845 

0.028423  (n  =  7.036n') 

3 

0.6721845 

0.0683353  (n  *  2.93n’) 

4 

0.791012 

0.126167  (n  =  1.585n') 

5 

0.8388 

0.173449  (n  =  1.1531n') 

6 

0.8531 

0.193578  (n  =  1.033n') 

7 

0.856412 

0.198812  (n  =  1.006n’) 

Recall  that  n  (the  replacement  rate  in  the  finite  spares 
case)  is  0.20.  The  repair  rate  n  *  required  to  achieve  the 
same  level  of  system  availability  as  achieved  in  the  finite 
spares  case  is  always  less  than  n,  but  must  converge  to  n  as 
n  gets  large.  The  ratio  between  n  and  n '  is  shown  in  the  last 
column  of  Table  2.2. 
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D.  DISTRIBUTION  OF  DOWN  TIME  DUE  TO  LACK  OF  PARTS 


In  a  mission  of  a  fixed  duration,  the  contribution  of 
the  nth  spare  part  cannot  be  determined  solely  by  looking  at 
the  point  availability.  This  is  because,  even  with  generous 
spares  support,  the  system  will  alternate  between  up  and  down 
states  as  the  system  fails  and  is  replaced.  Ono  indicator  of 
the  contribution  of  the  nth  spare  part  is  the  decrease  in  down¬ 
time  that  results  from  the  inclusion  of  the  spare  part.  In 
this  section  we  address  this  problem  by  deriving  the  distribu¬ 
tion  of  downtime  due  to  the  lack  of  spare  parts. 

Suppose  we  have  n-1  spares  available  for  a  component  for 
a  mission  of  duration  [0,t]  .  j^s  before  let  be  i»i.d.  random 

variables  (exponential)  representing  "up"  times  and  {R^}?^ 
be  i.i.d.  random  variables  representing  replacement  times. 

Let 


X 

1 


I 


i=l 


T 

i 


+ 


n-1 


l 

i=l 


R 

i  ‘ 


is  the  random  variable  representing  the  duration  from  time 
0  to  the  point  in  time  at  which  the  component  fails  for  the 
nth  time. 

If  an  nth  spare  (n+1  parts)  is  added  to  support  this  com¬ 
ponent  then  the  nth  spare  will  begin  to  function  after  a  re¬ 
placement  time  of  length  Rn.  Let  X2  =  X1  +  Rn'  then  the  actual 
contribution  of  the  nth  spare  to  the  availability  is  made 
after  time  X^. 


J 
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Let  Y  =  max( 0 , t-X2 ) .  Then  Y  is  a  random  variable  repre¬ 
senting  the  duration  of  "down"  time  due  to  the  lack  of  more 
than  n-1  spare  parts.  Let  k^  (t)  be  the  .pdf  of  the  random 
variable  X^.  Then 

k(h)(t)  =  [ (f*g) (n) ] (t) 

and  the  distribution  function  is  given  by 


Ly(x) 


P{  Y  <_x}  =  P{X2>_t-x} 


/  k(n) ( t ) d  i 
t-x 


For  example,  when  n  =  2, 


k(2)  (t) 


Uf*g)  (2)J 


(t) 


tf]2[(T  -f->e‘XT  +  (t  +|)e'nT] 


*<2)(x) 


dL  (x) 
dx 


k<2) (t-x) 


tfl 2  [  (t-x-  f )  e'1  +  (t-x  *  f)  e-”  1  *-*> ) 


for  0  <  x  <  t,  and 


i(2)(0)  =  /  k(2)(T)di  =  l|-l 2  [  /  (x-f)e'XTd 


t  “  t 

+  /  <x+f)e"nTdx] 


since  Y  *  0  when  X2  _>  t. 
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Thus,  the  p.d.f.  of  Y  when  n  =  2  is 


t  r9 1  2  r  i^"  ±  ^  2  x  “ x t  ,  # t  1  ,  2  x  ” n 1 1  y  —  r» 

lJ]  [(T  +  -T  I7)e  +  (n+T  +  ^)e  1  for  X  =  0 

A  n 


i{2)  (X)  =  J 


and 


W  [|]2[  <t-x -|)e‘X  (t_x)  +  (t-x +|)e~n  (t_x)  ]  for  0  <  x  <_t 


(2-7) 


Graphically, 


The  closed  form  expression  for  the  distribution  function 


of  X  for  arbitrary  n  is: 


l/n)(x) 


,8,nre-»(t-!!)  "-1 

[?i  ‘-TH^rn-  ii0 


I  (-D1 


(n-1) !  (t-x) 

(n-i-l)  ! 


n-i-1 


(-X) 


i+r 


n-l  _  ,,P.  -X (t-x)  n-r-1 

£  (-l)r  vn+r-1;  r.e 

Y=1  r !  6r 


e  '  ' r"  T  */  ,,j  (n-r-1) 

(n-r-1)  1 [  (  1]  "(n-r-j-i) 


(t.r..x,Ln  f  j  _ ] )  +  (— 1 ) n  t^-1  n{  — — t— 

(-X) j+1  5 


-n(t-x)  n-l  . 

-ITT  l,  <-D  • 


(n-l)  1  .(t-x)1*'1-1  nrx  (n+r-l)Pr  e~n(t~x) 
* (n-i-l) 2  (-x ) i+i  r-l  r !  6r  ' 


•  ( 


n-r-1 


I 

j-0 


(-1) 


i  (n-r-1)  !  . 
•■(n-r-j-U'j 


(t-x)^’^1 

. -11 


(2-8) 


If  we  define  (x)  to  be  1  -  l/n'  (x)  ,  then  the  expected  down 

time  in  the  interval  (0,t]  when  there  are  n  spares  is  given 
by 

/  (x) dx  , 

0 


and  the  expected  amount  of  down  time  prevented  by  adding  the 
nth  spare  part  is 


/fc  (L(n-1) (x)  -  L(n) (x) )dx  . 
0 
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E.  OPERATIONAL  AVAILABILITY  FOR  A  COMPONENT  WITH  SPARES  AND 

A  REPAIR  CAPABILITY 

In  the  previous  work  we  have  considered  the  case  in  which 
the  single-component  system  had  finitely  many  spares  and  no 
repair  capability.  Thus,  when  all  the  spares  are  used  a 
failure  will  cause  the  system  to  go  down  and  remain  down 
until  the  end  of  the  mission  deployment  period  at  which  time 
repairs  will  be  made  or  additional  spares  provided.  Because 
of  manpower,  equipment,  space,  and  training  problems  little 
actual  repair  is  done  during  say,  a  deployment  cruise  of  a  ship. 
Thus,  the  finite  spares/no  repair  problem  is  operationally 
realistic.  If,  however,  a  repair  capability  does  exist  during 
the  mission  deployment  period,  the  number  of  spares  required 
to  provide  a  specified  level  of  operational  availability  can 
be  reduced  substantially .  Figure  2.1  depicts  the  type  of 
situation  that  exists  when  spares  and  a  repair  capability  bot.. 
exist. 


Figure  2.1:  Failure/Replacement/Repair  Process 
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Several  researchers  have  worked  on  the  problem  in  which 
a  repair  capability  exists.  Barlow  and  Proschan  [  5]  consider 
the  case  in  which  the  system  consists  of  a  single  component 
with  one  spare,  instantaneous  replacement,  and  a  single  repair 
server.  They  obtain  the  mean  system  downtime,  the  Laplace 
transform  of  the  distribution  of  the  time  to  the  first  failure 
and  the  limiting  value  for  system  availability  for  general 
failure  and  repair  time  distributions.  Gnedenko  [19]  and 
Subramanian  [40]  consider  the  case  of  an  N  component  series 
system  with  a  single  spare  and  a  single  server  repair.  They 
obtain  the  Laplace  transform  of  the  availability.  Subramanian 
and  Ravichandran  [42]  extend  these  results  to  an  arbitrary 
number  of  spares.  Kumagai  [23]  considers  a  single-component 
system  with  n  spares  and  a  single  server  repair.  He  obtains 
an  symptotic  formula  for  stationary  availability  for  the  case 
of  exponential  lifetimes  and  general  repair  times.  Gaver  [17] 
considers  the  same  case  and  obtains  the  Laplace  transform  of 
the  probability  density  of  the  time  to  system  failure,  the 
mean  time  to  system  failure,  and  the  average  long-run  fraction 
of  time  during  which  the  system  will  be  available.  All  of  the 
above  research  assumes  the  replacements  are  instantaneous  when¬ 
ever  spares  are  available.  With  the  repair  capability  there 
are  essentially  infinitely  many  replacements  possible  during 
an  interval  of  time.  In  our  problem,  there  are  only  finitely 
many  replacements  possible  during  a  given  mission  and  the 
replacement  times  are  not  negligible. 
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A  similar  problem  has  been  considered  by  other  researchers 
with  interest  being  in  measures  of  effectiveness  more  closely 
related  to  supply  system  performance.  Schrady  [37 1  *  Richards 
[34],  and  Sherbrooke  [39]  have  examined  repairable  item  inven¬ 
tory  models  which  assume  that  some  failures  will  be  lost  to 
the  system  as  unrepairable.  Figure  2.2  is  modified  as  shown 
in  Figure  2.1  to  include  losses  to  the  system  and  a  source  of 
external  supply.  The  efforts  of  these  authors  have  been  concerned 


with  determining  the  number  of  spares  and  replenishment  and 
repair  policies  so  as  to  minimize  down  time  during  a  period. 

Suppose  that  a  single-component  system  is  supported  by  k 
spares.  When  the  system  fails,  replacement  takes  place  immediately 
if  there  is  a  spare  available,  and  the  failed  component  enters 
the  repair  process  immediately.  (Observe  that  this  is  equiva¬ 
lent  to  the  assumption  that  k+1  identical  units  are  in  parallel 
with  repair  of  failed  units.)  The  system  is  down  if  and  only 
if  the  number  of  units  in  the  repair  facility  is  k+1.  Let 
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q^(t)  be  the  probability  that  i  units  are  in  the  repair  facility 
at  time  t  for  i  =  l,2,...,k+l.  The  operational  availability 
is  then  given  by 


A(k) (t)  =  1  -  qk+1(t)  . 


For  the  special  case  k  =  0,  a  failure  causes  the  system  to 
go  down  and  remain  down  until  the  component  is  repaired.  The 
state  of  the  system  behaves  as  an  alternating  renewal  process. 
Analysis  is  mathematically  equivalent  to  our  initial  model 
with  infinitely  many  spares  and  non-zero  replacement  times. 

It  is  well  known,  (see  for  example  toss  [35])  that  the  opera¬ 
tional  availability  for  the  case  in  which  failure  times  are 
exponential  iid  with  parameter  X  and  replacement  times  are 
exponential  iid  with  parameter  n  is  given  by 

AM(t)  «  TTT  +  exp(-t(X+n)) 


On  taking  the  limit  as  t  -<-<»,  the  steady  state  operational 
availability  is 


lim  A(oo)  (t) 
t -*•'=» 


n  _  MTBF 

X+n  MTBF  +  MTTR 


This  is  "definition"  four  for  operational  availability  dis¬ 
cussed  in  the  Introduction.  We  see  that,  mathematically,  it  is 
correct  for  the  case  in  which  one  is  interested  in  the  limiting 
result  when  there  are  iid  exponential  failure  times,  iid 
exponential  replacement  times,  and  infinitely  many  spares  (or. 


41 


equivalently,  0  spares  and  a  repair  capability) .  For  k  >  0,  the 


determination  of  the  probability  q^+^(t)  is  much  more  compli¬ 
cated.  Choi  [8],  Ashar  [3],  Srinivasan  [40],  consider  this 
problem. 

In  the  next  chapter,  we  extend  our  analysis  to  multi- 
component  systems  with  finite  spares  support.  We  return 


to  the  situation  in  which  no  repair  capability  exists. 
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III.  POINT  AVAILABILITY  OF  TWO-COMPONENT  SERIES  SYSTEM 
WHEN  TWO  COMPONENTS  ARE  INTERDEPENDENT 

A.  INTRODUCTION  OF  SCENARIOS 

A  two-component  series  system  is  maintained  by  replacing 
failed  components  when  they  fail.  The  times  to  failure  for 
type  i  components  are  exponentially  distributed  with  parameter 
x ^  and  the  replacement  times  are  exponentially  distributed 
with  parameter  The  component  lifetimes  are  independent. 

With  respect  to  what  happens  to  the  surviving  component 
while  a  failed  component  is  being  replaced  the  following  five 
scenarios  are  defined. 

Scenario  1:  The  surviving  component  continues  in  service, 
and  if  it  fails,  its  replacement  begins  immediately  and  pro¬ 
ceeds  independently  and  concurrently  with  the  replacement  of 
the  other  failed  component.  All  failed  components  resume 
operation  as  soon  as  they  are  replaced. 

Scenario  2:  The  surviving  component  is  shut  down  until 
replacement  of  the  failed  component  is  accomplished. 

Scenario  3:  A  type-2  component  continues  in  operation  while 
a  type-1  component  is  being  replaced,  but  a  type-1  component 
shuts  down  while  a  type- 2  component  is  being  replaced. 

Scenario  4:  The  failure  rate  of  the  surviving  component  changes 
to  e^X^  ( 9 ^  <  1)  while  the  failed  unit  is  being  replaced. 

Scenario  5:  The  failure  rate  of  the  surviving  component 
changes  to  9^X^  (9^  >  1)  while  the  failed  unit  is  being  replaced. 


43 


Very  little  previous  work  has  been  done  for  the  problem 
in  which  the  system  components  are  not  independent.  Liang  [24] 
considered  2-component  series  systems  like  those  described  as 


scenarios  1,  2,  and  3  with  infinitely  many  spare  parts.  He 
obtained  the  limiting  availability  for  the  case  in  which  life¬ 
times  and  repair  times  are  exponential.  We  are  interested 
in  determining  the  point  availability  of  the  systems  for  both 
the  finite-spares  and  the  infinite-spares  cases. 

B.  SYSTEM  POINT  AVAILABILITY  WHEN  COMPONENTS  ARE  INDEPENDENT 

(SCENARIO  1) 

Let  A^(t)  be  the  availability  of  the  system  composed  of 
two  components  in  series  operating  under  scenario  1.  We  assume 
that  both  components  start  at  time  0  from  "up"  states.  Since 
the  system  is  a  series  system,  it  will  be  operational  at  time 
t  if  and  only  if  both  components  are  operational.  Since  the 
components  operate  independently,  the  probability  that  both 
are  operational  at  t  is  the  product  of  the  probabilities  of 
components  one  and  two  being  operational  at  time  t.  Let  a^(t) 
and  a2(t)  be  the  component  availabilities  at  time  t,  then 

A^t)  =  a1(t)a2(t) 

In  fact,  in  the  independent  case,  this  result  can  be  extended 
to  a  series  system  composed  of  any  finite  number  of  components. 
The  system  availability  of  a  k-component  series  system  (when 
the  components  operate  independently)  is  given  by 
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k 

A.  (t)  =  n  a.  (t) 

1  i=l  1 

This  result  can  be  extended  further  to  include  parallel  systems 
and  mixed  series/parallel  systems  provided  the  components  oper¬ 
ate  independently.  The  result  for  an  n-component  parallel 
system  is 

n 

(1  -  a.  (t))  =  n  (1  -  a.  (t)) 

i=l  1 

(The  unavailabilities  satisfy  a  product  rule.)  The  reader, 
familiar  with  mathematical  reliability  will  recognize  these 
results  as  identical  to  those  used  to  calculate  the  reliability 
of  a  series/parallel  system  configuration.  In  fact,  the  availa¬ 
bility  of  a  mixed  system  configuration  is  calculated  just  as 
the  system  reliability  is  determined.  The  component  point 
availabilities  simply  replace  the  component  reliabilities. 

For  the  case  in  which  components  1  and  2  in  series  both 
have  infinitely  many  spares  and  exponential  iid  lifetimes  with 
parameters  A ^  and  a2,  respectively  and  exponential  iid  replace¬ 
ment  times  with  parameters  and  n2,  the  system  availability 
is: 


Ax(t) 


(J!i_ 

y\1  +  n 


1 


+ 


*1 

xl  +  nl 


-(A.+n, )  t 
e  1  A  ) 


•  ( 


n2 

X2  +  n2 


+ 


X2 

A2  +  n2 


- (A_+n~) t 
e  *  *  ) 
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For  the  case  in  which  the  components  have  finitely  many 
spares,  expression  2-1  from  the  previous  chapter  should  be 
used  for  a^(t). 

C.  SYSTEM  POINT  AVAILABILITY  FOR  SCENARIO  2  (INFINITELY 
MANY  SPARES) 

Let  A  be  the  state  in  which  both  components  are  up; 

B  be  the  stae  in  which  both  components  are  down; 

be  the  state  tht  component  i  is  down  and  the  other 
component  is  up. 

Let  T^  be  the  lifetime  of  the  first  unit  of  component  1 
and  T^  the  lifetime  of  the  first  unit  of  component  2.  Let 
R^  be  the  replacement  time  of  component  i.  Then  the  time  to 
the  first  failure  of  the  system  is  given  by  min[T1»T2]  which 
has  pdf 

-(x  +x  )  t 

h(t)  *  (A1  +  x2)e 

The  recurrence  time  of  state  A  depends  on  which  component  fails 
first.  If  component  i  fails  first,  the  recurrence  time  is 
min[T1,T2]  +  R^ . 

The  conditional  pdf  of  recurrence  time  given  component  i 
failed  first  is  then  given  by 

t  -(A.+X,)s  -n. (t-s) 

r. (t)  =  /  (x.+x_)e  n.e  ds 

x  0 

The  probability  that  the  first  failure  is  of  type  1  is 
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P [X  <  Y] 


Similarly, 


/  P [X  <  Y |X  =  t] f (t)dt 
0 


I  P  (Y  >  t)f  (t)dt 
0 


»  -x_t  -x.t 
/  e  2  x.e  1  dt 

A 


X  1  +  X  2 


P  [ Y  <  X]  = 


X  X  • 

12 


Therefore,  the  pdf  of  the  recurrence  time  of  state  A  is 
given  by 


x,  t 


r  (t)  = 


L  Ul+X2)e 


1  "2  0 


-(X1+X2)s  -n, (t-s) 

n^e  ds 


x  2  t 


+  t — +y~  /  (x ,  +  x _) e 
A1  a2  0  1  2 


-(X.+X-Js  -n0(t-s) 

1  n2e  2  ds 


A  set  of  those  first  recurrence  times  forms  a  renewal 
process.  Let  h(t)  be  the  renewal  density  function  of  such  a 
process,  then  we  have 

t  __ 

A2(t)  =  FiF2  +  J  h(u)du.F1(t-u)F2(t-u)  .  (3-3) 
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Taking  Laplace  Stieltjes  transformations 

*  1  ■*  ■^‘1  n  n  1 

r  ( s )  =  x  ^  •  -  -  x 

1  s+X^+X  2  s+n^ 

n  X  -  +X  -  T1  - 

r  (s)  =  ■L  £  .  _ 

2'  '  s+X X+X 2  s+n2 
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s+n^  s+r^  s+a^+a2 


h* (s)  = 


r*  (s) 


h*(s)  = 


X1n^(s+n2>  +X2ti2(s+n^) 
(s+X^+x2)  (s+n^)  (s+n2) 


(s+X1+x2)  (s+r^)  (s+n2) 


(s+X^+X2)  fs+n^)  -  X^n-^  ( s-Hri 2 )  X2^2  ^s+ni^ 


xlnl^s+r|2^  +  ^2n2  ^s+nl^ 


(s+X^+X2)  (s+n^)  (s+r^)-*!^  (s+n2) -^2n2  ^s+ni^ 


il{  (h*FxF2)  {t) }  =  h*  (s)  • 


s+x  ^+x  2 


where 


*lnlS  +  ^2n2S  +  ^Xl+X 2^ nln2 


s(s+X-^+X2)  [s  +(X^+X2+n^+n2)  s+tn-j^+X-j^+Hj^) 


Xlnl  +  X2n2 


( x 2.+x 2 )  n^n. 


(s+x1+x2)  (s-r (s-r2)  s  (s+X1+X2)  (s-r.^) 


“  (xi+x2+ni+n2^  ±  r^Al+x  2+nl+n2^  ”4  (n^ri2+x1n2+n^X  2> 


(k  =  1,2) 


since 


/*“!  ^  ^ 2 n2  ^l+*2^  nln2 

^  {s+x^+x 2>  (s-r^)  (s-r2)  +  sTi+T^TrpTi-r^TTs-T^T 

n,  n,  a ..  n ,  +x _ n _-n  .  n-  -(A  +x  )t 

"  ^  +  (A^a^Hx,^^  e 


^Xlnl  + X2n2^rl  +  ^Xl+X2^ nln2  rlfc 
+  r^(A^+A2+r^)  (r^-rj)  e 

r2^Xlr'l  +  X2n2^  +  ^Xl  +  X2^nln2  r2fc 
+  r2(Xi+X2+r2)  <r2_rl)  6 


So 


A2(t) 


-(Xx+X2)t 


+  nin2 
rlr2 


X lnl + X 2n2 ~ nln2  -(Aj+* 

+  ^Xl+X2Hrl^  ^Xl+X2+r2^  S 


( x  1  n  i+x  2^  2  ( X  ^+X  2 )  ri 2  r^t 

+  r1(x1+x2-fr1)  (rx-r2)  6 

r2  (Xini+X2ri2^  +  ^Xl+X2^ nln2  r2t 
+  r2Ul+X2+r2)  (r2"rl)  6 

It  can  be  shown  that 


Xlnl  +  X2n2  ”  nln2 

(x1+x2+r^) (xx+x2+r 


-1  . 


Finally 
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+ 


A2  (t) 


nln2 

rlr2 


(X1n1  +  X2n2*ri  +  ^xl  +  X2^nln2 


rl (rl  +  +  x2) (ri 


r2) 


( X 1  ^  1  +  X2n2^Y2  +  ^X1  +  X2^nln2 
r2(r2  +  x1  +  x2)  (r2-rx) 


rlfc 


.r2fc 


(3-4) 


D.  SYSTEM  POINT  AVAILABILITY  FOR  SCENARIO  3  (INFINITELY 
MANY  SPARES) 

Let  A(t)  be  the  state  of  the  process  that  both  components  are 
in  service  at  t,  B(t)  be  that  both  components  are  out  of 
service  at  t,  and  D^(t)  be  that  component  #i  is  out  of  service 
with  the  other  component  in  service. 

Setting  up  a  set  of  Kolmogorov's  forward  differential 
equations , 


A’  (t) 

DJ(t) 
d2  ( t) 

B'  (t) 


=  ”  ( X  1+x  2 )  A  ( t)  +  n1D1(t)  +  n2D2(t) 

=  XxA(t)  -  (n1+X2)D1(t)  +  n2B(t) 

=  X2A(t)  -  i2D2(t)  +  nxB(t) 

-  x2Dl(t)  “  (n1+n2)B(t) 


Taking  Laplace-Stielt jes  transformations  yields 


SAMs)  -  A  ( 0 )  =  -(Xj+X2)  A*  (s)  +  niD*(s)  +  n2D2*(s) 
sD*  (s)  -  D^O)  =  A1A*(s)  -  (n1  +  n2)DJ  (s)  +  n2BMs) 
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.  • — —*• 


ps*!r 


X2A*(s)  -  n2D*(s)  +  niB*(s) 


sD*(s)  -  D2(0)  = 

sB*(s)  -  B  (0 )  =  A2D*(s)  -  (n;L  +  n2)B*(s) 


using  initial  conditions  and  solving  a  system  of  linear 
equations  for  A*(s),  we  have 


A*  (s) 


nln2 


*l*2nl 

P-i  ->  +  X  ..  n  *5  +  n  ^  ^  \  TT 

12  12  12  ni+n-.+A, 


'1  '2 


1 

s 


2  (  n  2+X  2-r>]_ ) 

(n2+^2^  Ui  +  nl~*2_n2^ -^iX2^ 


s+r. 


3  2  2  2  2 

-r2+  (2n-j_+2n2+^  2^  r2"*  ^nln2+nl+n2+nl*  2+n2'V2^r2+nlr|2;'l2+nlr|2+' 

-r2(r1-r2)  <r3-r2> 


(s+r2J 


3  2  2  2  2 

~r^+  (2n^+2". 2+x  2)  r^~  (3n^r^+r^+n 2+ni^2+n2^2^ r 3+nln2*  2+n ln  2+ 

-r3 (ri“r3) ^2~r3^ 


(s+r3) 


where 


rl  =  (n2+^2) 


53 


and 


I  2  2 

U^+.^+Sn^+r^)  ± k  (xi+X2+2nl+n2^  ~4  ( n]_A  i+nix  2+n2X  l+nln2+nl  ^ 


(k  =  2,3),  are  the  eigenvalues  of  the  infinitesimal  transition 
matrix . 

Finally,  the  closed  form  solution  would  be 


A3(t)  = 


nln2 


xlX2nl 

n  !  n  2+X  !  n  2+ n  !  X  2'tT]T^+r^ 


n  ix  2  ( ri2+X  2~nl^ 


- (X  2+n2^ 


(x2+n2)  (n-[_n2+x j  X2+nlX2~nlXl~nl  ^ 


3  2  2  2  2  2 

-r2+ (2n1+2n2+X2^  r2~  ^nln 2+ni+n2+nlX  2+n2X 2^  r2+  ^ nln2X 2  +  nln2+nln2^ 

'(r“-r  )  (r3~r2)  ~ 


* 

> 

i 


3  2  2  2  2  2 

-r3+(2n^  +  2n2+X2^r3~^riin  2  +  ni+n2+n  1X  2+n  2X  2^  r3+  ^  nln2X  2+nln  2+nln2^ 


(-r3) (r1-r3) (r2-r3) 


.-r3fc 


(3-5) 
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Notice  that  all  eigenvalues  of  the  infinitesimal  transition 


matrix 


D  = 


<X1+X2> 


(X2+nl) 


- (xl+n2) 


( n-^+n  2^ 


rl ,  r^  and  r^  are  all  real.  To  see  this,  consider  the  square 


1 2  auu  an  real.  to  see  uus,  consiaej 

root  expression  found  in  equations  for  rj  and  r-j. 


( xi+x  2+^ni+n2^  ^  nlxl+nlx2+n2Xl+r|l  n2+r,l^ 


—  +n2  ^ri2+^X2~^Xl)  =  /^Xl_ri2^  +^2  ^xi+X2+^n2^ 

which  is  real  since  all  and  take  positive  values.  Also 
it  is  easily  seen  that  all  eigenvalues  are  positive.  As  a  re¬ 
sult  of  this,  the  availability  is  known  to  decay  exponentially 
with  time . 

E.  LIMITING  RESULTS  FOR  THE  AVAILABILITIES  FOR  SCENARIO  1,  2, 
AND  3 

It  is  simple  to  find  limiting  availabilities  (lim  A. (t)  =  A. ) 

by  applying  results  from  the  limiting  behavior  of  Markov  process 

[10]  .  Let  tt  denote  limiting  probability  of  state  s  in  the 
s 

embedded  Markov  Chain. 

For  scenario  2,  the  transition  probability  matrix  of  the 
embedded  Markov  Chain  having  states  A,  ,  D2  is 
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Xl+X2 


0 

0 


Xl+X2 


From  the  identity  n  = 


ttP  ^ ,  i  •  e . 


IT  =  ’T_  +  7T_ 

A  Dx  D2 


D1  x i+x  2  A 


D_  A.+i  A  ' 

£  x  Z 

and  from  the  limiting  result  of 


A_  =  lim  A_ ( t)  = 


nAuA 


t-Va> 


2  Va^d/d^D^I^ 


where  u  ,  u  ,  and  un  represent  the  mean  sojourn  times  of 
a  Uj.  u2 

states  A,  D, ,  and  D_,  respectively.  We  have 


A  \ ^+A  2 


1  X1  1  x  2  1 
+  r - — — . - 7T»  +  t - — — - it. 


A  ^]_  +  ^2  Xl+X2  nl  ^  Xl+X2  n2  ^ 


Hi  n 


12 


nln2+Xln2+X2nl 


(3-6) 


which  we  also  obtain  from  eq.  (3-4) .  We  have  used  the  facts  that 


xl+x2 
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and 


For  scenario  3  the  transition  probability  matrix  of  the 
embedded  Markov  Chain  having  states  A,  D^,  D2,  B  is 


nl+X2 


Xl+X2 


nl+n2 


xl+x2 


nl+n2 


nl+X  2 


For  the  identity  of  n  =  ttP. 


- -  1T_  + 

nl+X2  ^1  ^2 


X1  +  n2 
A ,  +A  „  A  n, +n. 


12 

X. 


12 


B 


2  t  '1 

.  IT  -  1  TT  T 

x  i  ^x  ^  A  n  1  +n  o  B 


1TA2 


'1  "2 


T11+X2  ^ 


Solving  this  system  of  equations  in  terms  of  tta  gives 
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A ^  (n^+n2) (n1+A2> 
nD^  nl^Xl+X2^nl+n2+X2^  ^ 

A  2 ( A  ^+A  2+n i+n2^ 

*D2  =  (A1+A2)  (n1+n2+A2>  ’'A 

A^A2  (  nl+n2^ 

nB  =  n1(A1+A2) (ni+n2+A2)  *A 

and  using 


Xl+X2  ' 


nl+A 2  ' 


D. 


and 


v 


B 


1 

nl+n2 


the  limiting  result 


A 


3 


lim  A,(t) 

t-i-oo  J 


’Va 


1IAUA  +  "n  Hn  +  "r.  P 


D1  D1 


D2  d 


becomes  (by  substitution) 
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nln2(nl+n2+X2) 


nln2^nl+n2"t’X2^+Xlri2^nl+ri2^+rilX2  ^ X]_+^2+nl+ri2^  +n2XlX2 


nln2 


X1X  2nl 

nln2  +  xLn2  +  nxX2  + 


(3-7) 


Note  the  relationship  between  the  limiting  availability  expressions 
and  "mean  system  up  time"  and  "mean  system  down  time".  The 
limiting  availability  of  the  system  is 


A  = 


mean  system  up  time 


mean  system  up  time  +  mean  system  down  time 


(3-8) 


According  to  this  formula  (3-8) ,  expressions  (3-1) ,  (3-6)  and 
(3-7)  can  be  written  as 


Xl+X2 


Xl+X  2 


l  tXl  +  +  X1X2 

A. ^+A 2  n 2 


(3-9) 


X 1  +  X  2 


1  *  l  *2  1 

z  +  /_£  +  _ ±)  z 

x^+x2  n2 


(3-10) 


and 


A-  = 


Xl+X2 


1  ,  1  [Xl|nl+n2)  ,  A2Ul+X2~>'ril'fn2^|  X1X2 

X ^+x 2  x  +X  ri  ^  ( n n2+X 2 )  n2  (n^tr)2tx2^  ^"tri 2"tx 2) 

(3-11) 
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respectively.  We  can  see  from  (3-9) ,  (3-10)  and  (3-11) 
that  mean  system  down  time  under  scenario  1  is 


1  X1  ,  X2  ,  X1X2] 
*  1+* 2  ^i  n2  r|in2 


i 


that  of  scenario  2  is 


Xl+X2 


and  that  of  scenario  3  is 


1  Xl(nl'fn2) 
A^+A2n^(n^+n2+A2) 


X2  ^nl+n2+Xl+X2^ 


X1X2 


(n1+n2+A 


Let  Y i  be  a  random  variable  representing  system  down  time 
under  scenario  i,  then  from  the  above  observations  we  can  see 
the  relationship 


st  st 


Y1  i 


Y3  i  Y2 


(3-12) 


which  is  equivalent  to 


EtYj,]  >  E[Y3)  _>  E[Y2] 


Also  from  (3-9),  (3-10)  and  (3-11)  we  can  see  the  relationship 


Ax  1  A3  <  A2 


(3-13) 
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To  contrast  these  series  system  results  with  those  of 


parallel  system  consider  the  limiting  unavailability 


B  =  lim  B  ( t) 


UB  _  ^B^B _ 

PBB  "  rBV(lD;D1+1D2PD2+VA) 


where  wBB  denotes  mean  recurrence  time  of  state  B. 
Unavailability  in  this  parallel  system  is 


_ mean  system  down  time _ 

mean  system  down  time  +  mean  system  up  time 


1 

nl+n2 


1 

nl  +  n2 


+ 


1 

nl+n2 


nln2, 

X1X2 


(3-14) 


From  (3-14)  we  can  see  that  "mean  system  up  time"  in  the  two- 
component  parallel  system  is 


1  rnl  nln21 

ril+r,2  X1  X2  X1X2 

which  is  very  large  compared  with  those  in  the  series  systems. 

Notice  that  a  symmetry  exists  in  the  formulae  (3-9)  and 
(3-14)  . 

F.  POINT  AVAILABILITIES  WHEN  COMPONENTS  SHARE  A  COMMON  POOL 
OF  SPARES 

Under  the  assumption  that  the  two  series  components  share 
a  common  pool  of  finitely  many  spares,  we  can  still  use  the 
Markov  process  results. 
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Define 


q^j/lws)  =  probability {Xn+1  =  k,tn+1  -  tR  =  s/Xn  =  j} 
with  states  A,  D^,  D^,  B  for  scenario  i,  then 


0 

fF 

fF 

0  ' 

0 

fF 

fF 

gF 

0 

0 

fG 

fG 

'  q2  = 

q 

0 

0 

gF 

0 

0 

*  q 

0 

0 

0 

gG 

gG 

0  . 

'  0 

fF 

fF 

0  - 

0 

fF 

fF 

0  ' 

gF 

0 

0 

fG 

gF 

0 

0 

fG 

9 

0 

0 

0 

'  qk  = 

gF 

0 

0 

fG 

.  0 

gG 

gG 

0  , 

0 

gG 

gG 

o. 

k  =  4,5 

*  —  0  A  t 

where  f(t)  =  eXe  .  In  scenario  2,  state  B  does  not  exist. 

Each  element  of  the  above  matrices  represents  the  p.d.f. 
of  sojourn  time  of  state  j  multiplied  by  the  probability  that 
the  next  state  the  process  will  visit  is  k.  For  example,  in 
the  element  fG(s)  (=  Xe  (x+n)s)  could  be  written  as 
(n+x)e  ( n+A  )  s  .  _A_  •  This  means  that  at  time  s  the  process 
leaves  state  and  the  probability  that  state  B  will  occur 

before  state  A  occurs  is  t-t— . 

X+n 

If  the  system  is  available  at  t  with  the  kth  spare  alive  at  t 
then  2k  transitions  have  been  made  before  t  since  each  spare  used 
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requires  2  transitions.  Hence  availability  of  the  system  under 
scenario  i  with  n  common  spares  is 


A|n) (t) 


A|n-1) (t) 


t 

I  IqJ 

r\ 


(2n) 


1.1(8) 


■Fz (t-s)ds 


=  A|n_1)  (t)  +  {  Iq^210^ ^*F2}  (t)  (3-15) 

where  [qf^n^).  1  is  (1.1)  element  of  2n-fold  convolution  of  the 
matrix  q^.  This  recursive  formula  can  also  be  written  as 


A{n)  (t)  =  F2(t)  +  .IiAq^2j)]i.1(s)-F2(t-s)ds 

*  F2(t)  +  [  (tq{2j))1  ^F2}(t)  (3-16) 

j3*  J- 

Notice  that  { [qf2-^ ] ^  ^*F2}^  rePresents  the  marginal  contribu¬ 
tion  to  the  availability  of  the  jth  spare.  For  example,  if 
n  -  2 

A{2)(t)  =  F2  ( t)  +  2[fF*gF*F2]  (t)  +  4[  (fF)(2)*(gF)(2)*F2J  (t) 

+  4  IfF*fG*gF*gG*F2]  (t) 

A<2) (t)  =  F2 ( t )  +  2[fF*g*F2]  (t)  +  4[ (fF)  (2)*g(2)*F2]  (t) 
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A^2) (t)  =  F2 ( t)  +  [fF*gF*F2] (t)  +  [fF*g*F2] {fc) 

+  [(fP) (2)*g(2)*F2] (fc)  +  2[(fF) (2)*g*gF*F2] (t) 

+  [  (fP)  (2)*  (gP)  (2) *F2] (t)  +  [fF*fG*gG*gF*F2] (fc) 

+  [fF*fG*gG*g*F2] (t) 
or  alternatively 

a{2)  (t)  =  A^15  (t)  +  4[(fF)  U)*(gF)  (2)*F2]  (t) 

+  4  [fF*fG*gF*gG*F2] 

A^2)  (t)  -  A^15  (t)  +  4[(fF)  (2)*g(2)*F2]  (t) 

A^3) (t)  =  A^X) (t)  +  [ (fP) (2) *g(2)*F2] (t) 

+  [ (fF) (2) *g*gF*F2] (fc)  +  [ (fF (2) * (gF (2) *F (2) ]  (t) 
+  [fF*fG*gG*gF*F2] (fc)  +  [fF*fG*gG*g*F2] {fc) 

Here 

(a{2)  (t)  -  a{15  (t)]  =  {[qf4)]L1*F2}(t) 

is  the  marginal  contribution  to  the  availability  of  the  second 
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Fig.  3.1  Sample  paths  for 


h.l  F  }(t) 


Closer  inspection  of  the  term  (t)  is  useful. 

In  order  for  the  system  to  be  alive  at  t  with  the  second  spare 
operational  the  sequences  of  states  the  sample  paths  for 
this  term  will  follow  can  be  categorized  as  following  eight 
cases.  Since  a  symmetry  exists  by  the  change  of  role  between 
component  #1  and  #2,  the  illustration  of  four  cases  is  enough 


Probf  Case  II)  =  Prob{Case  III}  =  [fF*fG*gF*gG*F2]  (fc) 


ProbfCase  IV}  =  [ { f F ) (2) * (gF) (2) *F2] 


(t) 


Notice  that  by  the  summetry 


2  *Prob{ Case  II  U  Case  IllUCase  IV} 

=  2[(fF) (2)*(gF) (2)*F2] (t) 

+  4 [fF*fG*gF*gG*F2] 

-  [ (f*g*F) (t) ] 2 

which  represents  the  probability  that  starting  from  state  A 
at  time  0  the  process  will  be  found  in  state  A  at  time  t 
with  exactly  one  spare  replaced  in  each  component  some  time 
before  t. 

To  see  the  last  equivalence,  recall  from  (2-3)  that 


{ [f*g*F] {t)}' 


r  9  -nt  .  9  .  -At  9  -At. 2 
[—7  e  +—  te  - 7  e  ] 

5 


Now  computing  the  other  terms  gives 


[ (fF)  (2)*(gF)  (2)*F2] 


(t) 


=  £-x[- 


( n  A ) 


(S+2A)  (s+n+A) 


-] 


(nA) 


(s+2A ) 3 (s+n+A ) 2 


A  , 

(  S  +  2  A  )  + 


BCD 
- - J  +  -(s  +  n+")  + 


( S  +  2  A )  (s  +  2A ) 


( S+n+A ) 
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where 


A  = 


Inverting  the  Laplace  transforms 


KfF)  (2)*(gF)  (2)*f2]  {fc) 


3e  -2xt 

?  e 


29' 
6  : 


te‘2Xt  + 


02  t2 

7*^ 


-2X 


e~ (n+x ) t 


e2  .  -(n+x)t 


Similarly 


[fF*fG*gF*gG*F2] (fc)  =  £1( 


_ ( n  x )  _  , 

(s+2X ) 2 (s+2n) (s+n+x)2 


5e 2  -2x t  02  .  -2xt  92  -2nt 

—  —  e  +  — t  te  +  — r  e 


26' 


45 


46 


.  e2  -(n+x)t  .  e2  -(n+x)t 
+  —r  e  +  — r  te 

6  6 J 


By  substitution  we  have  the  above  equivalence. 

G.  EXTENSIONS  OF  SCENARIOS  FOR  MORE  THAN  TWO-COMPONENT 
SERIES  SYSTEMS  WHEN  SPARES  ARE  SHARED 

Since  all  the  components  in  this  system  share  the  same 
spare  parts  we  can  formulate  a  semi-Markov  kernel  [9]  over 
the  state  space  according  to  the  given  scenario. 
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Taking  an  example  of  a  three- component  system  will  illus¬ 
trate  the  extension.  Let  D. ..  be  the  state  that  component  i, 
3  and  k  are  not  in  service,  and  consider  the  following  exten¬ 
sions  of  scenarios. 

Case  1:  Extension  of  scenario  2,  i.e.,  if  any  of  the 
three  components  fails  then  the  others  will  be  shut  down. 

VJe  have  a  semi-Markov-kernel  q2  having  states  A,  D^,  and 


*2 


g 

g 

g 


Case  2j  Extensions  of  scenario  3. 

Consider  the  following  four  subcases: 

Case  2-a:  If  component  #1  fails  first  the  others  will  be 
shut  down,  but  the  operations  of  component  #1  are  not  affected 
by  failures  of  others,  then  the  semi -Markov  kernel  q2a  will  be 


0 

fF2 

fF2 

fF2 
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0 

0 

0 

g 

0 

0 

0 

0 

0 

0 

0 

gF2 

0 

0 

0 

fFG 

0 

fFG 

0 

gF2 

0 

0 

0 

0 

fFG 

fFG 

0 

0 

gFG 

gFG 

0 

0 

0 

0 

fG2 

0 

gFG 

0 

gFG 

0 

0 

0 

fG2 

0 

0 

gFG 

gFG 

0 

0 

0 

fG2 

0 

0 

0 

0 

gG2 

gG2 

gG2 

0 
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over  the  states  A,  D 


1' 


D2 '  D3 ' 


12'  13' 


D23  and  D123 


Case  2-b:  If  component  #1  or  component  #2  fails,  component 
#3  will  be  shut  down,  but  a  failure  of  component  #3  does  not 
affect  the  operations  of  component  #1  or  component  #2.  In 
this  case  the  semi-Markov  kernel  is 


f  0 

fF2 

t-h 

N) 

fF2 

0 

0 

0 

0 

gF 

0 

0 

0 

fG 

0 

0 

0 

gF 
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0 

0 

fG 

0 

0 

0 

q3b 
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0 

0 

0 

0 

fFG 

fFG 

0 

0 

gFG 

0 

gFG 

0 

0 

0 

fG2 

0 

0 

gFG 

gFG 

0 

0 

0 

fG2 

,  o 

0 

0 

0 

„F2 

gG 

gG2 

gG2 

o  , 

over  the  states 

A,  Dlf 

°2 ' 

°3 '  °12 '  °13' 

°23 

and  D223* 

Case  3-< 

c:  The  component 

#3  shuts  down  only 

when 

both  of 

components  #1, 

and  #2 

are  not  in 

service 

.  The  semi 

-Markov 

kernel  in 

this 

case  is 

I 

f 

0 
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0 

0 

0 
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0 
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over  the  same  states  as  in  case  3-b. 


Case  3-d 

Is  If 

any 

two  of 

the  components  fail 

the 

other 

component 

shuts 

down 

.  The 

semi- 

Markov 

kernel 

would 

then 

0 

fF2 

fF2 

fF2 

0 

0 

0 

0 

gF2 

0 

0 

0 

fFG 

fFG 

0 

0 

gF2 

0 

0 

0 

fFG 

0 

fFG 

0 

gF2 

0 

0 

0 

0 

fFG 

fFG 

0 

q3d 

0 

gG 

gG 

0 

0 

0 

0 

0 

0 

gG 

0 

gG 

0 

0 

0 

0 

0 

0 

gG 

gG 

0 

0 

0 

0 

.  0 

0 

0 

0 

gG 

gG 

gG 

0 

over  the  same  states  as  in  case  3-b. 

Once  we  set  up  a  semi-Markov  kernel,  A^11'  (t)  is  obtained 
by  eq.  (3-16) .  Actual  evaluations  get  harder  as  n  increases 
in  eq.  (3-16) .  In  cases  such  as  these  approximation  formulae 
or  availability  bounds  would  be  necessary. 
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IV.  COMPUTATIONAL  FORMULAE  FOR  THE  AVAILABILITIES 
OF  DIFFERENT  SYSTEM  CONFIGURATIONS 


A.  INDEPENDENT  TWO-COMPONENT  SERIES  SYSTEM 
1 .  Case  of  No  Sharing  of  Spares 


{A2,n2} 

Denote  by  A^n,m^  (t)  the  availability  of  the  system 

when  component  #1  has  n  spares  and  component  #2  has  m  spares. 

Both  components  start  their  operations  from  "up"  states. 

Denote  by  A^  (t)  {c  =  I,II)  the  availability  of 
c 

component  c  when  j  spares  are  available.  Then 

AUi,m)  (t)  ,  AU)  (fc)  .  A£n)  (fc)  (4-1) 


spare)  . 

I  would 

Agn+l'm)  (t)  =  {A^n)  (t)  +  [  (f1*g1)  (n+1)*F1]  (t) 

=  A^n'm)  (t)  +  t(f!*gi)  (n+1)*F1]  (t) -A^}  (t) 

.4-2) 


Suppose  we  add  one  more  spare  to  component  I  (n+lth 
The  marginal  increase  in  the  availability  of  component 
be  [  (f^*g^)  ^n+^*F^]  (t)  •  Hence 
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Also, 


A(n'm+1) (t) 
s 


Ain'm)(t)  +  [(f2*g2) (m+1)*F2] (t) -A^n) (t) 


(4-3) 


When  we  add  (1.1)  spares  to  (n,m) , 


A^n+l,m+l)  (t)  =  A^n,m)  (t)  +  [  (f1*g1)  (n+1)  *FX)  (t)  *A^}  (t) 

+  t(f2*g2) {m+1)*F2] (t) • A^n) (t) 

♦  [ (*!**!>  In+1,*F1] (t) • [ (f2*g2)  ""+1,*f2I (t) 

(4-4) 

The  second  term  of  the  R.H.S.  of  eq.  (4-4)  is  the  contribution 
made  by  the  (n+1)  st  spare  of  component  I  to  the  system 
availability i  the  third  term  is  the  contribution  of  component  II;  and 
the  last  term  is  made  jointly  by  the  (n+1)  stand  tm+1)  st  spares. 

2 .  Case  of  Sharing  of  Spares 


1  11 

1  I  1 

{A,n)  { x,n) 


Assume  there  are  n  spares  available  for  both  components. 
Hence,  if  component  I  uses  k  (0  <_  k  <_  n)  spares  during  [0,t] 
then  component  II  can  use  up  to  (n-k)  spares.  Let  A^^  (t) 
be  the  probability  that  the  system  is  available  at  time  t  with 
(i,j)th  spares  alive  at  t.  Then 
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A<l)(t) 


a'°-®U> 


-2Xt 


A(0'0}(t)  +A{1,0}(t)  +  A{0,1} 
s  s  s 


=  A^u; (t)  +  2[f*g*Fj (t) -F(t) 


(t) 


A{0,0}(t)  +A{1,0}(t)  +A{0,1}(t)  +A{0,2}(t) 
s  s  s  s 

+  A{1,1} (t)  +  A{ 2,0} (t) 
s  s 


a11)  (t)  +  2[(f*g)  (2)*F]  ...  *F(t)  +  ( [f*g*F]  ..  . }: 


One  schematic  expression  for  the  possible  events  which  comprise 
A (t)  is  shown  below  for  n  =  4  -(n+1  =  5)  . 


Component  1:0,  1,  2,  3,  4,  5* 


Component  II  :f  0 


Therefore, 


Ai5)(t)  =  Ai4>(t)  +  A*.°’5}(t)  +Al1,4}(t)  +  a12,3} 


+  A^3,2}(t)  +  A{4,1}(t)  +  A^5,0}(t) 
s  s  s 
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(t)  +  2F(t> • [ (f*g) 


*i5’  <*> 


=  A 


(4) 


(5>*p 


F] 


(t) 


+  2 [  f *g*F )  (t)  *  t  (f*g)  l4)*F]  (t) 


+  2 [ (f  *g) (2)*F] (t) -I (f*g)  (3)*F] 


(t) 


A<4) (t)  +  IQ[(f*g)  lx,*F] (t) -t(f*g) 


(5_l)*F] 


(t) 


where 


[  (f*g)  (0)*F]  (fc)  =  F(t)  . 


In  general. 


A(n+i) (t) 


=  Afn) (t )+  l  C (f*g) '*/ *F] 
s  1=0 


(n+1) 


(t) 


•  [  (f*g)  (n+1'l)*F]  (4-5) 


B .  INDEPENDENT  TWO-COMPONENT  PARALLEL  SYSTEM 
1 .  Case  of  No  Sharing  of  Spares 


1 

1 

1  ( A  .  „n. 

L} 

1  . 

11 

1 _ 

[ 

{ A2  r n2  ^ 

Denote  by  Apn'm^(t)  the  availability  of  this  parallel 
system  when  (n,m)  spares  are  available.  Then 
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AUi,m)  (t)  =  A(n)  (t)  +  Ata)  (fc)  _  A_(n)  (t)  .A(m)  {fc)  (4-6) 

If  we  add  one  more  spare  to  component  I ,  then 

AUi+l,m)(t)  =  {A(n)  (t)  +  [(f1*g]_)  (n+i)*F]  (t)> 

+  A*™3  (t)  -  { A^n)  (t) 

+  t(f1*g1)  (n+1)*F13  (t)}-A^3  (t) 

A^n+l,m)(t)  =  AU,m)  {t)  +  [  (*1*g1>  (n+1)  *^1  (t) 

-  [(f1*g1)  (n+1)*F1]  (fc) -a£J3  (t)  (4-7) 

A^n,m+1) (t)  =  A(n,m) (t)  +  [ (f2*g2) (m+1) *f2J (t) 

-  [(f2*g2) (m+1)*F2] (t) -A^n) (t) •  (4-8) 

The  marginal  contribution  made  by  the  (n+l)th  spare  of  component 
I  is  A^n+3-,m3  (t)  -  A^n,m3  (t)  ,  and  that  of  the  (m+l)th  spare 

*r  C 

of  component  II  is  A(n'm+1) (t)  -  A(n'm3  (t) .  When  we  add 

P  P 

(1.1)  spares  to  (n,m)  spares, 

AOi+l,m+l)  (t)  =  A^n,m)(t)  +  [  (f L*gL) (n+1)*F1] (t) 

■n-A^ttn  *  [  <fj*g2)  <”+1)*r2!  (t| 

• (1  -  A^P)  (t) ]  -  [ (f1*g1)  (n+1)*F1]  (t) 
•{(f2*g2llm+1,*r2](t)  .  (4-9) 
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» 


2 .  Case  of  Sharing  of  Spares 


U,n) 


ii 


U,n) 

Denote  by  (t)  the  availability  of  this  parallel 
system  when  n  spares  are  available  for  both  components.  As 
in  IV. A. 2, 


Ap0,‘t>  - 

=  2e-Xt  -  e'2Xt 

II 

rH 

+  A{1'°}(t)  +  A{0'^ 

P  P 

(t) 

= 

4J 

r*-> 

o 

o 

+  2A{1,0} (t) 

P 

= 

A<0) (t)  + 

2{{f*g*F]  {t)  +  F  ( t)  - 

(f *g*F] 

In  general 

A(n+D  (t) 

P 

=  A^n) (t) 

+  (  l  5  A  (t) 

i=0  ? 

(4-10) 

where 


Aptl,n+l"i}  (t)  =  [  (f*g)  (i)*F]  (t)  +  [  (f*g)  (n+I-i)*F]  (t) 


-  [ (f*g)  (l)*F]  (fc)  • [ (f  *g)  (n+l-l) *F] 


(t) 


(4-11) 
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and 


[  (f*g)  (0)*F]  (t)  =  F ( t) 


or 


A^n+1) (t) 


A(n) (t)  +  2  [  A 

P  i=0  P 


(?) 


{ i,n+l-i } 


(t) 


if  n  is  even 


A(n) (t)  +  2  l  A 

P  i-0  p 


(2§i) 


{ i ,n+l-i } 


(t) 


, n+1  n+l. 
t— 5— /—j— J 

+  A  1  1  (t) 


if  n  is  odd 


where 

A^j,k}(t)  =  [  (f*g)  (j)  *F]  (fc)  +  [  (f*g)  (k)*F]  ((;) 

-  [  (f *g)  (j)*F]  (t) • [ (f*g)  (k) *F ]  (t)  . 

C.  INDEPENDENT  MIXED  CONFIGURATION 


U3 » n3 ) 

[Block  I]  [Block  II] 
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(4-12) 


(4-13) 


;#V 


Let  (t)  be  the  availability  of  Block  I  which  is  one- 
component  subsystem  consisting  of  component  #1. 

Let  (t)  be  the  availability  of  Block  II  which  is 

two-component  parallel  subsystem  consisting  of  components  #2 
and  #  3 . 

Let  be  the  availability  of  this  mixed  configuration 

with  (n,m,£)  spares  are  available.  Then, 

A^n'm'i)( t)  =  A^n)  (t)  • 

■  Ajjj  (t)  •  [A j 2]  (t)  +  A[3]  (t)  “  A[2]  (t) 

•  A | (t) 1  (4-14) 

where 

(k) 

Aji] (t)  =  availability  at  t  of  component  #i  (i  = 1,2,3). 

When  we  add  one  more  spare  to  one  of  those  components, 

we  have 

A^n+l,m,a) (t)  =  A(n+l)(t,.A(m,t)(t) 

=  (A^n)  (t)  +  t(f1*g1)  (n+1)*F1l  (t)  }.aJ“'M  (t) 

=  AjJn*m*fc)  (t)  +  [  (f1*g1)  (n+1)*F1l  (t) 

•  A^,l)  (t)  (4-15) 
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^n.m+1,  i)  (fc)  =  A^)(t).  Ai^.O(t)  +  [  (f2*y2)  (m+1)*F.J 


2  (t) 


U  "  A^j  (t)  ] 


*  A^n,m'°(t)  +  A^n)  (t)  •  [  (f  2*g2)  (m+1)  *F2]  (fc) 


11  -  A | ^ | (t) ] 


(4-16) 


^n.nwt  +  l)  (t)  =  A^n,m'°(t)  +  A^n)  (t)  •  [  (f,*g,)  "  '  *'  *F 


( ^  +1 ) *p  i 

3  y3;  F3] (t) 


[1  -  A 1 2 1  (t)  ] 


(4-17) 


Suppose  component  #2,  and  component  #3  are  identical 
such  that  they  share  p  spare  parts.  Then 


<'pl  it) 


aJ"’ (t) 


A^ht) 


(4-18) 


Adding  one  more  spare  to  the  first  component 


A^n  +  1 #  P )  ( t )  = 


A^n)  (t)  +  [  (fj/g^  (n+1)*F1]  (t)  } -A^P)  (t) 


=  A^n‘p)(t)  +  ((f1*g1)(n+1)*F1](t).A;[P)(t) 


(4-19) 


and  adding  one  more  spare  to  the  second  block  gives 
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4n'p+1)(t) 


(p+1) 


A^n-p)  (t)  +  A^n)  (t)  • [' *1  ' ( [f2*g2)  {1) *p2]  (t) 
+  !<f2*52>(P+1'i)*fr21(t) 

[  (f2*g2)  *  ^  *^2^  (t)  ■  ^  tf  2*g2^P  ^*F2^(t)* 


(4-20) 


The  marginal  contribution  made  by  the  (p+l)th  spare  of  block 
II  is 


A<n'p+1)(t) 


(t)  . 


D.  DEPENDENT  TWO-COMPONENT  SERIES  SYSTEM  FOR  THE  "SYMMETRIC 
SHUT  DOWN"  SCENARIO 

1 .  Case  of  No  Sharing  of  Spares 


{X1,n1>  {A2,n2} 


Component  #1  starts  its  operation  at  time  0  with  n 
spares  and  component  #2  begins  with  m  spares.  The  system  is  avail¬ 
able  if  (1)  the  original  units  of  components  1  and  2  survive  to 
time  t  or  (2)  the  (k, i. )  th  spares  (k  =  1 ,  .  .  .  ,  n,  i  -  1,  .  .  .  ,m)  are  alive 
at  time  t.  In  the  case  of  (2)  the  following  situations  can  occur. 
Suppose  that  the  first  failure  occurs  at  time  s  <  t  and  the  unit 
is  replaced  at  time  u  (s  <  u  <_  t)  by  its  first  spare  part  letting 
the  process  start  all  over  again  from  u  with  one  fewer  spares. 
Figure  4.1  shows  this  process  for  the  case  in  which  the  first 
failure  occurred  in  component  1. 
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So  the  following  recursive  formula  follows. 

A^n,m)  (t)  =  F1F2(t)  +  nf1F2*g1)*A^n“1,m)  ]  (t) 

+  [  (f  2fl<*2)  *^{2n'm~1)  1  (t)  •  (4-21) 

Computational  formulae  can  be  derived  by  looking  at 
the  process  differently.  Let  A^'^ft)  be  the  availability 
under  scenario  2  with  (i,j)th  spare  alive  at  t.  Looking  at 
the  first  couple  of  terms. 
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A^3,0)(t)  =  A<2'0)(t)  +  A^  3 ' 0 }  ( t) 


(2,0) 

’2 


(3)*=r  F 

1*2  yly  FlF2J(t) 


=  A;-'*"  +  [(f-F  *g.)  '  ' *F. F, ] 


a^1,;l)  (t) 


A^1'0*  (t)  +A2°'1}(t)  +  A^1 * 1 ^ ( t) 


A2  *  ^  +  2Fl*g2^  *F1F2^  (t) 


+  (x) [ (f1F2*g1) * (f2Fl*g2} *F1F21 (t) 


we  have  (^)  since  either  component  can  fail  first. 


A^2,1) (t) 


A^1,X)(t)  +  A22,0;(t)  +  A22,1}(t) 


=  A 


(1,1) 


(t)  +  (  (f1F2*g1)  (2)*F1F,] 


1  2  (t) 


+  (l11<flF2*n)(2)*<f2Fl*92)*FlF2'(t) 


A^2,1) (t)  +A23,0>(t)  +A23,1)(t) 


A2  r  ^  +  ^flF2*gl*  *F1F2*  (t) 


+  (]_)  £  (fiF2*g1)  (  (f2Fl*g2}  *F1F21  (t) 


In  general. 
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A^n,m) (t) 


A<n'm‘1)(t)  +  [  (“t1) 

z  i=0 


[  (flF2*gl) 


(i) 


(f2Fl*g2) 


(m) *F. F_ ] 


r  2  (t) 


(4-22) 


or 


Ain,m)  (t)  _  A^n-l,m)  (fc) 

m  .  .  ... 

+  jl0  <Y}  •nfiF2*gi)(n)Mf2Vg2)(3)*FiF2](t)  (4-23) 

where 


[h(g) *F) 


(•) 


=  F 


(•) 


Either  of  eq.  (4-22)  or  eq .  (4-23)  is  preferable  to  eq.  (4-21) 
for  computational  purposes. 

To  evaluate  eq.  (4-22)  or  eq.  (4-23)  we  need  to 
evaluate  the  term  [  ( f  ^  *  (f  2F1*92^^  *F1F2^  (t)  * 


illl'iV’l1  ll,*(f2r1*g2)  (t)> 


<3) 


(X1n1)  i  (Aj7^  ^ 

(s+^+A2)1+^+^(s+n^)^(s+n2)g 


For  (i,j)  =  (1,1)  we  have 


(4-24) 


(s+x 

Let  a  = 

Thus , 

l(f, 


Then 


(X^^n^)  (X2n2^ 


2^+X 2 )  (s+n^)  (s+n2> 


^  1 1  ^  (X2n2}  [  (s+X1+\2) 


B  C 

+  - 2  +  - 


( S  +  X  ^  +  X  2 )  (S  +  X1  +  X2) 


(s+n^)  Ts+n^T 


—  X^  —  X2#  b  n2"*^2_™^2*  then 


A  = 


2  2 

az  +  b  +  ab 

373 

a  b 


B  =  - 


a  +  b 

. 377 

a  b 


C  ab 


D  =  - 


3  ' 

a  { n  2“  ’i  i ) 


E  = 


b  (n2“n1) 


F2.g1)M£2F1*g2)*F1F2)(t) 


, ,  ,  ,,  x  r 3  +b  +ab  12 

(X1n1) (x2n2^  ^  373  6 

5L  D 

a+b  ’"(X1+X2)t  1  .t2  -(xi+x2)t 

-  te  +  IE*“  e 

a  b 

-  -TT^ - 7  e  +  -  e  ]  • 


a  ( n  2~ n  ^ ^ 


b  {r\2~r\^) 


As  another  example  consider  the  case  (i.»j)  =  (2,1) 
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-(X  +X  )t  -(x.+x  )t 
=  (X1n1)z(X2n2HAe  x  +  Bte  A 

2  -(A  +A-)t  3  -(A.+A  )t 

+  C-^~  e  A  *  +  D-Jy  e 

-n-it  -n.t  -n~t 

+  E-e  x  +  Fte  +  Ge  } 


where 


A 


(4b+3a)(  ( 3b+2a)  (2b+a) -3ab)-ab{  5  (2b+a)  +3  (3b+2a) -3  (a+b)  } 


B 


C 


D 


E 


F 


1  (3b+2a) (2b+a)-3ab 

2*  ,4h3 

a  b 


2b+a 

aV 


a  b 


4  ( n 2  —  n  -^ )  -a 

TFT  r? 

a  'r'2"ni 


a 


a  ( n  2“ ^ i ) 


and  a  -  ' i~A 2 '  b  =  n2~^l_X2' 
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As  (i,j)  gets  larger  it  becomes  more  difficult  to 
invert  the  Laplace  transform.  In  these  cases  it  is  often 
reasonable  to  make  some  assumptions  which  simplify  the 
expressions . 

In  many  real  situations  the  replacement  rates 
are  very  large  compared  to  the  failure  rates  x . .  When  this 
h-appens  the  expressions  involving  e  are  dominated  by  the 

-Xit 

expressions  involving  e  and  good  approximations  for  the 

availability  can  be  obtained  by  ignoring  the  terms  involving 
“nit 

e  .  The  formulae  can  be  simplified  significantly  without 
-nit 

setting  e  to  zero  if  a  single  replacement  rate  nj_  =  n2  ~  n 
can  be  assumed.  The  expression  for  availability  using  the 
common  replacement  rate  n  is  much  more  tractable. 

It  can  be  shown  that  using  max{n^>  for  the  common 
replacement  rate  will  give  an  upper  bound  of  the  true  availa¬ 
bility  and  min(n^}  a  lower  bound.  A  weighted  replacement  rate 
may  be  used  for  a  good  approximation  too. 

In  any  case,  let  the  common  replacement  rate  be  n,  then 
eq.  (4-24)  becomes 


iU£1F2*g)  <ll*(£2F1*g)  '^‘FjP.,]  (t) ) 


(4-25) 


Taking  the  inverse  transformation  gives 
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-(X1+X2)t 


UfjF^g)  (i)*(f2F1*g)  (:i)*f1f2]  (t) 


,lX2r'1+j  ,  ti+j  .  nr  (i+j+r-l)Pr  t1+j_r 

■^r+3  *n+jrr  rIi  _1  r j  5r  TiTFFH 


.  ,  v  i+j+1  xix2nl  jrt1+j'1  1+?_1  (i+j+k)  Pk 

.i+j-k-l  . 

t  _ _ ,  -nt 

.+i-k-l)!1 


(i+j-k-l) 


where 


o  = 


n-X  x~X  2 


and 


nPk 


=  n (n-1) . . . (r-k+1)  = 


n! 


(n-k)  ! 


2 .  Case  of  Sharing  of  Spares 


{X,n>  { A,n) 


Using  the  same  notations  as  in  IV. A. 2, 


A^X) (t)  =  A20,0}(t)  +  2A^1,0}(t) 


=  A^0)  (t)  +  2  [fF*g*F2] 


(t) 


(4-26) 
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A<2>(t) 


A^1} (t)  +  2A^2,0}(t)  +  A21,1}(t) 
A^1}  (t)  +  (2  +  2)  [ (fF*g)  (2) *F2] 


A<3)(t) 


A^2/ (t)  +  2A^3'0}(t)  +  2A^2,1}(t) 


and  in  general 


A{22)  (t)  +  [2  +  2(3)]  .[(fF*g)  (3)*f2]  } 


A^n) (t) 


A^n~1)  (t)  +  2n[ (fF*g)  (n) *F2] 


(t) 


(4-27) 


which  can  also  be  obtained  from  (3-15) .  And 


[ (fF*g) (n}*F2] 


(t) 


a  n  j.n  n  ,  p  .  n-r 

[«’  Ji(  11  ~ Ti”' 


n-1 


.  ,>n+l ,  e11  ft0-1 

+  {  X)  [7HTT]  [TivrTTT 

0 


n-1  .  ...  P.  .  n-£-l  , 

.  v  (n+s.)  k  t  ,_-nt 

i  l  (n-£-1 ) I Je 

l-l  il 6  in  11  • 


(4-28) 


where 

e  =  r\\  , 

6  =  n  -  2x  . 


E,  A  NUMERICAL  EXAMPLE 

Consider  a  two-component  series  system.  We  assume  for 
simplicity  that  the  components  are  identical.  We  consider 
two  cases:  1)  spares  are  shared  by  both  components,  and 


2X  t 


90 


I 


2)  each  component  has  a  separate  spares  pool  with  no  sharing 
allowed.  With  the  pa-ameters  used  in  the  previous  example, 
x  =  yjj-,  n  =  and  t  =  90,  we  obtain  the  results  shown  in 
Table  4-1  from  equations  (4-5),  (4-28),  (2-2),  (2-3),  and 
(4-26)  . 


Table  4-1:  Computation  Results 


Spares  Shared 

Spares  not  Shared 

n  |  Scenario  1 

Scenario  2 

(nl,n2i 

Scenario  1 

Scenario  2 

Scenario  3* 

0  i  0.00248 

0.00248 

(0,0) 

0.00248 

0.0024788 

0.00248 

1  0.01914 

0.02293 

(1,0) 

0.01081 

0.012704 

0.0.081 

2  '  0.07154 

0.09404 

(0,1) 

0.01081 

0.012404 

0.01270 

,  3 

0.17750 

0.23416 

(1,1) 

0.04713 

0.05848 

0.05786 

4 

0.31973 

0.41220 

(2,1) 

0.10036 

0.1288 

5  |  0.46724 

0.56961 

(1,2) 

0.10036 

0.1288 

0.58768 

0.67178 

(2,2) 

0.21371 

0.2659 

7 

0.66670 

0.72204 

(3,2) 

0.31074 

0.4208 

_ _ 

(2,3) 

0.31074 

0.4208 

1 

(3,3) 

0.45180 

0.5668  | 

The  following  computations  show  how  the  result  for  A 


(1,1) 

3 


(t) 


was  obtained. 


A<°'0}(t) 


=  F2(t) 


A^'0}(t) 


[f *g*F ] (t) *F (t) 
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A^ ' 1  (t)  =  [fF*g*F2](t) 

A^lfll;t)  =  2[  (fF)  (2)*g*gF*F2]  (t)  +  [ f F* fG*gF*gG*F2 ]  ( t) 

+  [f F*fG*gG*g*F2 ] (t) 

and 

2  [  (fF)  k2)  *g*gF*F2]  (fc)  =•  2  ( n  \  )  2  t  A  +  Bt  +  |  t2  )  e  2Xt 


where 


A  = 


(2n-3\r  -  (  n-\ )  (n-2\) 
( n-2\ ) 3 ( n- X ) 3 


B  =  - 


(  2  n  -  3  \  ) 
1 - 


( n-2  X )  ( n- \ ) 


T  ' 


c  = 


(n-2  \ )  (  n-  \ )  ' 


D  = 


3  ' 

(n-\)X 


E  =  - 


( n-2  \ )  X 


[fF*fG*gF*gG*F2]  (t)  =  (n\)2[(- 


5  .  t  \ „“2\ t 

T  + 


4  ( 3-  \ )  2  ( n-  \ ) 


T  + 

(n-\)  {  n  —  X  > 


4  <  n- X ) 
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and 


[fF*fG*gG*g*F2] (fc) 


=  <n\)*[(- 


4  <n-A)  3(n-2A)  2 ( n-A ) 2 ( n-2A ) 


)e 


-2 


1  - ( n+\ ) t  ^  1 

e  + - 


(v-03\ 


e~ nfc  +  — i 


n\  ( n  —  2 \ ) 


4n  ( n-A ) 


e“2nt] 


Thus  we  have 


A^Ut) 


A^0,0f(t)  +  A3il,0}(t)  +  A3{0,1^(t)  +  A^1,1}(t) 


=  0.00248  +  0.00833  +  0.01022  +  0.03683 


0.05786  . 
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V.  COMPARISONS  OF  POINT  AVAILABILITIES  FOR  DIFFERENT  SCENARIOS 

In  this  chapter  we  continue  to  consider  the  point  availa¬ 
bilities  of  two-component  series  systems  which  operate  in 
accordance  with  scenarios  one  through  five.  Because  of  the 
complications  that  arise  when  trying  to  obtain  explicit 
mathematical  formulae  for  the  point  availabilities  in  some  of 
the  cases  examined  earlier,  we  attempt  in  this  chapter  to 
obtain  order  relationships  among  the  point  availabilities  of 
the  five  scenarios.  Such  order  relationships  allow  bounds  to 
be  placed  on  the  point  availabilities  for  the  scenarios  that 
generate  the  most  computational  problems.  Intuitively,  we 
would  expect  that  the  scenario  point  availabilities  would 
satisfy : 


Ajn,m)  (t)  ^  A_(n 'm)  (t)  ^  A^n  'iTl)  (t) 

since  the  surviving  component  remains  in  operation  in  scenario 
one  (and  is  therefore  subject  to  failure)  while  a  failed  com¬ 
ponent  is  being  replaced;  the  surviving  component  remains  in 
operation  sometimes  in  scenario  three  (depending  on  which  com¬ 
ponent  failed) ;  and  the  surviving  component  never  remains  in 
service  in  scenario  two.  Intuitively,  the  more  the  components 
are  in  service,  the  more  they  should  fail  and  the  lower  would 
be  the  system  point  availability.  Similarly,  one  would  expect 
tha  t 
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(n,m) 

1 


(t) 


A|n,ra) (t) 


(n,m) 

2 


(t) 


and 


A<n'm)(t)  <  A<n'm)(t) 

This  chapter  attempts  to  prove  these  relationships. 

A.  COMPARISONS  BETT'TEEN  SCENARIOS  1,  2,  3 

In  chis  section  we  consider  comparisons  among  the  point 
availabilities  for  two-component  series  systems  described  by 
scenarios  one,  two,  and  three. 

Let  us  first  develop  a  few  results  that  will  be  useful 
in  deriving  the  order  relationships  among  the  point  availa¬ 
bilities.  Consider  a  one-component  system  with  n  spares 
available.  Let  A^n' (t)  be  the  availability  of  this  system 
given  that  the  system  started  up  at  time  0  with  n  spares. 

Let  (t)  be  the  availability  at  time  t  given  that  the  system 

started  down  at  time  0  with  n  spares  available.  The  following 
corollary  proves  that  A(t;n)  is  at  least  as  large  as  B(t;n) 
for  all  t-0  and  n.  That  is,  the  probability  that  the  sys¬ 
tem  will  be  found  in  an  up  state  at  time  t,  is  greater  when¬ 
ever  the  system  begins  in  an  up  state  than  if  the  system 
begins  in  a  down  state. 


0  to 


Fig.  5.1  IVo  Processes  for  A'n^  (t)  and  P'“ 
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Corollary  5-1 


A(n) (t)  >  B(n  (t)  t  >  0,  n  =  0,1,2,... 


Proof  s 

Let  s  be  the  time  of  the  first  failure  of  the  component. 
The  system  will  be  operational  at  time  t  provided  either  s  >  t 
or  s  <  t  and  the  system  is  resurrected  and  operational  at  time 
t.  In  the  latter  case,  at  time  s,  the  probability  of  being 
up  at  t  is  (t-s)  . 

A(n) (t)  *  F(t)  +  f (s)B(n) (t-s)ds  (5-1) 

0 

Similarly,  for  B^(t)  the  first  transition  (replacement)  must 
occur  at  some  time  s  <  t  at  which  point  the  probability  of 
being  up  at  time  t  is  A^n-1^  (t-s).  Thus, 

B (n) ( t)  -  /t  g(s)A(n_1) (t-s)ds  .  (5-2) 

0 

In  order  to  show  that  A^  (t)  21  B^  (t)  we  re-express  (5-1) 
and  (5-2)  as  follows: 

A(n)  (t)  =  F(t)  +  /  fG(s)B(n>  (t-s)ds  +  /  gF(s)A(n>  (t-s)ds 

0  0 

t  __ 

-  /  gF (s)F (t-s)ds  (5-3) 

0 


and 
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(t-s)ds 


B(n)(t)  «  /  fG(s)B(n) (t-s)ds  +  /  gF(s)A(n,-1) 

0  0 

Equations  (5-3)  and  (5-4)  were  determined  heuristically . 
Verification  of  the  equations  is  provided  by  showing  (5-3) 
to  be  equivalent  to  (5-1)  and  (5-4)  to  be  equivalent  to  (5-2) . 
We  do  this  by  showing  the  Laplace  transforms  to  be  identical. 
By  the  uniqueness  theorem  of  Laplace  transforms  this  proves 
that  the  expressions  are  equivalent. 

The  Laplace  transform  of  (5-3)  is  given  by 


A<“>*(s) 


1 

s+x 


s+x+n 


B<*>*(s) 


s+x+n 


A<n>*(s)  - 


s+X  +n 


1 

s+x 


which  simplifies  to  be 


A<n)* 


(s) 


X 

S+X 


*  [3) 


9 


which  is  the  Laplace  transform  of  (5-1) . 

Similarly,  the  Laplace  transform  of  (5-4)  is  given  by 


B  (n)  * (s) 


X  n(n)» 
s+x+n 


(s) 


+ 


n  A (n-1) * 
s+x+n 


(s) 


n 

s+X 


A (n-1) * 


(s) 


which  is  the  Laplace  transform  of  (5-2) .  Thus  (5-1)  is  equiva¬ 
lent  to  (5-3)  and  (5-2)  is  equivalent  to  (5-4) . 
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V  > 


Now,  in  (5-3) 


F  (t)  -  /  gF  (s)  F  (t-s)ds  *  e~(n+X)t  >  Q 
0 

and  from  (2-1) 

A(n>  (t)  >.  A(n_1)  (t)  for  t  >  0 

Thus, 

A(n)  (t)  =  /fc  fG(s)B(n)  (t-s)  +  J*  gF(s)A(n)  (t-s)ds  +  e"(X+n)t 

0  0 

>_  /  fG(s)B(n)  (t-s)  +  /  gF(s)A(n-1)  (t-s) 


Now  consider  a  two-component  series  system  operating 
under  different  operating  scenarios. 


{ X  ^ ^ }  { X  2  #  n2  J 

Denote  by  A^n,m^ (t)  the  availability  of  this  system  operating 
under  scenario  i  when  the  first  component  starts  to  operate 
from  the  "up"  state  with  n  spares  available  and  the  second 
component  starts  to  operate  from  the  "up”  state  with  m  spares 
available.  Then  the  following  theorem  holds. 


Theorem  5.2. 


A|n,m)  (t)  ^  A;(n,m)  (t)  <  A^n'm)  (t) 


t  >  0,  n,m  =  0,1, 2, 3,  ... 


The  proof  of  Theorem  5-2  is  inductive.  First  consider 
(n,m)  *  (0,0)  and  note  that 

A<°'°>(t)  -  A<°'0)<t)  -  A*0'0’*)  -  e”  (A^+Aj)t 


Next  note  that  for  (n,m)  =  (1,0) 


A^'0)  (t) 


[Fx(t)  +  (f1*g1*^1)  (t)]F2(t) 

P72(t)  +  l*l*2**lT2*ir2'  (t) 
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p^tt)  +  i«1F2*g1*F1r2]  (t) 


A?'0'*) 


FiF2(t)  +  lfip'>*giF9*?r1F?l 


1  2  yl  2  12*  (t) 


since 


F  ( • )  <  1 


A^tt) 


A^1,0) (t)  <  A*1,0) (t) 


Likewise,  we  have  for  (n,m)  *  (0,1)  the  relationships 


A^0#1) (t)  =  F1F2(t)  +  If 2Fl*g2Fl*FlF21 (t) 


A^0,1)(t)  =  F1F2(t)  +  [f2F1*g2*F1F2] 


(t) 


A^0'1) (t) 


*  FlF2(t)  +  tf2Fl*92*flF2] (t) 


and 


A^0'15  (t)  <  A^0,1)(t)  = 


A^'^U) 


Assume  that 


A^'^tt)  <  A^i,j) (3)  <  A^'j)  (t)  ,  t> 


and  for  all  (i,j)  such  that  i  <  n,  j  <  m,  then  we  will  show  that 


A^n,m)  <  A<n'm> (t)  <  A^n,m) (t)  ,  t  >  0  . 
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Before  continuing  with  the  proof  of  Theorem  5-2  we  define 
the  following  random  variables: 

Tj  :  Time  to  failure  of  the  first  unit  of  component  j 

Rj  :  Replacement  time  of  the  first  unit  of  component  j 

X(t)  :  State  of  the  process  at  t  (0  or  1) 

Vj(t):  State  of  the  component  j  at  time  t  (0  or  1) 

/V  o  \ 

We  also  let  '  '  be  the  availability  at  time  t  starting  from 
state  ;  (i  =  1,2)  with  (k,fc)  spares  available,  and  let 
T  =  minfT^T^,  T^  +  R^  “  S^,  and  T2  +  R2  =  S2.  Furthermore,  we 
define  the  following  seven  cases: 

Case  1:  T  =  T^  *  u^,  0  <  u^<  s^  <_  t  and  component  2  has  k 

failures  in  (u^s^l  with  the  kth  spare  alive  at  s^ 
So  at  s^  the  system  is  in  state  A  with  (n-l,m-k) 
spares  available  (Fig.  5-2) . 

Case  2:  T  ■  T^  -  u^,  0  <  u^  <  s^  <_  t  and  component  2  has 

(k+l)  failures  in  (u^,s^l  with  the  (k+l)th  replace¬ 
ment  underway  at  s^.  So  at  s^  the  system  is  in 
state  D2  with  (n-l,m-k)  spares  available  (Fig.  5-3). 

Case  3:  T  *  T^  »  u^O  <  u^  <  s^  <_  t  and  component  2  has  no 
failures  in  (u^s^  (Fig.  5-4). 

Case  4,  5  and  6  are  the  symmetric  analysis  of  cases  1,  2 

and  3,  respectively  obtained  by  changing  the  roles 
of  the  two  components  (Figs.  5-5,  5-6,  5-7). 

Let  pL(s,k)  be  the  probability  of  the  event  taking  place 
over  the  interval  (0,s^l  described  in  case  1.  Similarly,  let 
P2 (s^jk) ,  p3(s1),  p4(s2,k),  p5(s2,t),  and  P6(*2)  be  the 
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probabilities  of  the  events  taking  place  in  (0,s^l  described 
in  cases  2  through  6.  From  the  figures,  we  see  that: 


Pl(Sl,k)  *  j  [f1F2(u){g1[(f2*g2)  (}°*F2]  } 


u-Q  1  2  1  2  '2,i(«1-u)' 


p2(»i,k)  =  /u=i0I£1f2(u)(91[f2*g2l(,t,*£2*52n(Si.1 


p3(sl>  *  “iV®!1  (.1) 


p4<s2-1>  -  /  nlf2Fl(“)lg2I(fl*9l»U)*?l>>(,,-u)d'1 


Ps<32'1)  *  /a.0If2rl(“,(«2"fl*»l>“)**l*Sl>>(.2-«)a» 


P6(S2)  =  [f2Fl*g2] (s2) 

For  scenario  1,  cases  3  and  6  are  included  in  cases  2  and 
5;  therefore. 


A^n'm) (t)  =  F1(t)F2(t)  +  / 


t  m-1 


s1=0  k=0  A  A 

+  J1"  “l1  P2(s,,k) -A^-irm-k)  (t_s  )ds 

s1=0  k=0  2  1  iD2  11 

+  J  I  P4(s,,i) •A{n“£,m'1) (t-s,)ds. 

s2=0  *=0  4  2  1  24 

+  f  "‘I  Ps<«2»*>  •A{D"4'm"1)  (t-s,)ds. 

s2«0  z-0  3  2  iUl  2 


l  p1(s1,k)  .A|n_1'ro-k)  (t-s-^ds 


(5-5) 


likewise,  for  scenario  2,  cases  1,  2,  4  and  5  are  not  possible; 


therefore, 

.  (n,m) 


(t)  - 


F^tjF^t)  +  J  P3(s1)  •A^n’1,m)  (t-sx)dSl 


+  /q  P6(S2) *A2n'm'1) (t"S2)dS2  (5-6) 


For  scenario  3,  cases  4  and  5  are  not  possible  and  case  3 
is  contained  in  case  2.  We  obtain 


Ain,m)  (t)  =  F.  (t)F-(t)  +  /*  "'I1  Pl  (s.  ,k) -A^“  (t-ajds 

x  *  0  k=*0  x  1 


(n-l,m-k) 

3  \ 


+  /*  p  (S  ,k)  .A^_1'm"k)  (t-s.  )ds. 

0  k«0  i  1  J  2  1 


+  /  p6(s2) •A^n,m"1) (t-s2)ds2  (5-7) 


Let  us  now  compare  A^n,m) (t),  A2n'm) (t),  and  Ajn,m) (t)  using 
(5-5),  (5-6),  and  (5-7).  We  have. 
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<  f.F2+/tmi  P1(s1,k)A{n“1#m'k)(t-s1)ds] 
1  *  0  k=0 

0  k=Q  2 

+  f  “j1  P4(S  ,t)A<n-t'm-1)(t-3  )d. 

o  z*o  q  45  A 

♦  f*  "j1  P  (Sj.tlA}""1'”"11  (t-s  Ids. 
o  Jt-o  5  2  1  2 


<  F.F,  +  J*  mf  P  (s  ,k)A(n“1'm“k)  (t-s  )ds 
-  1  ^  0  k=0  -1  A  1 

♦  Z*  “j1  P  (»  ,k)A<g-1'B-lt,(t-.l)d,: 

0  k=0  2 

+  /fc  ^  P4(s2,UA}n'It‘1>  (t-s2)ds2 
0  l**0 

+  J*  V  P5(s2,UA{n'n,_1)  (t-s2)ds2 


P,P,  *  f  T  P1(s1.k)A<n-1'm-k) (t-s^ds 
x  *  0  k«=0  1  x 

+  /*  T  [;'srW<1','k',Mildl 

0  k=0  2 

♦  S„  (s2)  A"'”'1’  (t-Sl)dS 


<  F1F2  +  P1(s1,k)A|n"1'in"k)  (t-s2)ds2 


vs 


+  /fc  T  P2(S  ,k)A^"1,m'k)  (t-Sl)dS 

0  k=0  z  x  ±u2  ■L 

s2=sU  ^ 

=  Ajn'm) (t)  . 


1 


The  first  inequality  comes  from 

A^1#j)(t)  >  A^'j)(t) 

j 

which  is  immediate  result  of  corollary  5-1.  The  second  inequality 
comes  from 


A(i'j) 

A1 


(t) 


>  A 


(k,i) 

1 


(t) 


for  (i,j)  >  (k,  £)  and  t  >  0.  The  third  inequality  comes  from 
the  fact  that 

n~l  (  a  \  _  _ 

I  [  (f  *9i )  *  (F. +f i *Gi ) I t  \  =  P{Number  of  complete  cycles 

X,=0  1  1  1  it;  for  C0mp0nent  l  in  (0,t) 

1  n-1 }  ^  1 

for  n  =  0,1, 2, 3,...  .,  and  the  last  inequality  comes  from  the 
assumption  that 


A1 


(t) 


(t) 


for  (i,j)  <  (n,m)  and  x  ^  0.  It  is  easy  to  see  from  the  same 
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facts  used  above  that 


.  (n,m) 
A3 


(t) 


<  F.F  (t)  +  f f  F  *a 

1  r  1£  2  'T'1  irr  2  gl  A3  J  (t) 


+  [f2F1*g2A^n'm_1)] 


(t) 


<  FxF2  +  [f1F2*g1*A^n'1'm)] 


(t) 


+  [f2F1*g2*A^n'm_1)] 


(t) 


A<n'm) (t) 


This  completes  the  proof  of  Theorem  5-2. 

B.  COMPARISONS  BETWEEN  SCENARIOS  1,  4  AND  5 

The  following  corollaries  are  useful  in  comparing 
scenarios  1,  4  and  5. 

Corollary  5-3. 

Consider  two  single-component  systems.  Let  the  first 
have  exponential  iid  lifetimes  with  parameter  A^  and  exponen¬ 
tial  iid  replacement  time  with  parameter  n .  Similarly, 
let  the  second  have  exponential  iid  lifetimes  and  replace¬ 
ment  times  with  parameters  a2  and  n,  respectively.  If 
A 2  L  Aj then 

F x  (t)  F"2  ( t) 

F1(t)+[f1«5](t)  1  F2(t)+[f2*G] (t) 
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Proof : 


F^t) 


-X  .  t 
.  1 


lfi*sl  (t) 


X.  -x.t  . 
T“Ie  1  -e*"ntl 


n-x 


n-x. 


T~tl~e 


- (n-X  . )  t 


] 


Now,  let 


z  (t)  = 


n-x 


,  -(n-x,)t  x,  -(n-x.)t 

[1-e  1  ]  -  -4-[l-e  1  ] 


n-x. 


1 1  \  — (n— X  —  )t  x,  -  (n-x1 )  t 

-  Irrr-e  "  rrv-«  1  ] 


n-x2  n-x^  n-x2 


n-x. 


Taking  the  derivative  wrt  t,  we  see  that 


,  -(n-X,)  -(n-X.)t 

ft*(t>  "  V  “  le  1  ~  0 


since  X2  _>  x^. 

Since  1(0)  =  0  and  l( t)  is  increasing  in  t,  z  (t)  _>  0,  t  ^  0. 

Thus 


Fx(t) 


[fl*G](t) 


F2(t) 


[f2*Gl (t) 


and 

F,(t)  F,(t) 

- — - - -  _>  - = - — -  Q.E.D. 

F1(t)+[f1*G] (t)  F2(t)+[f2*G] (t) 

Note  that  F^(t)  +  [f^*G]  is  the  probability  that  the  first 
unit  of  the  ith  component  will  either  survive  to  time  t  or  will 
fail  in  (0,t]  and  not  be  replaced  by  time  t.  If  we  define  the 
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random  variable  N^(t)  to  be  the  total  number  of  cycles  (up 
time  plus  replacement  time)  of  component  i  in  (0,t],  then 
PtN^t)  =  0]  =  Fri ( t) -i- [ f i*G]  .  In  the  following  results  we 
use  this  interpretation  for  Corollary  5-3. 

Corollary  5-4. 

Suppose  all  the  assumptions  described  in  Corollary 
5-3  hold.  Let  (t)  and  (t)  be  the  availabilities 

of  the  first  and  second  system  respectively  and  assume 
also  at  time  0  that  both  systems  start  to  operate  from 
up  states.  Then, 

A<n)  (t)  >_  A^}  (t)  t  >  0  . 

Proof  of  Corollary  5-4 . 

We  have 

A<n)(t)  -  F^t)  +  JiHf1*g)(k)*F1J(t) 

and 

<«  ■  F2(t)  +  J0((f2*s)  <k>*F21  ,t) 

Let  X.^  be  the  random  variable  which  describes  the  sum  of 

uptime  and  a  replacement  time  for  component  i.  Since  x. 

st 

and  n1  =  n2»  we  have  X1  >_  X2  which  is  equivalent  to 


(5-8) 


(5-9) 

an 
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Prob{X1  >_  t}  >_  Prob{X2  >  t} 


or 


ProbtN^Et)  <_  k}  _>  Prob{N2(t)  <_  k}  (5-l( 


where 


Prob{N.  (t)  <  k}  -  l  [(f*g)  U)*(F.+f  *G) ]  . 

i  z  _q  i  x  i  ir; 


(5- 1_. 


Define 


(f^g) (kl (•> 


which  is  the  density  function  for  the  sum  of  k  cycles  for  com¬ 
ponent  i.  Rewriting  Eq.  (5-11)  yields 

k  t  .  .  F. (t-s) 

P(N.(t)  <  k}  =  l  j  hf  l)  (s).[P{N.(t-s)  »  0} - i - 

1*0  0  PU^tt-sJ-O} 

Ef,*G] . 

*  P(Ni(t-3)  -  0  )  ■p1M.(t-3i,,6)id3  <5‘12) 


Def ine 


and 


[fi*G](t-s) 
pTnT (t-s)  = b) 


*  Y 


i  dn 


(t,s) 


Then  we  have 


Yi  up  (t's)  +  Yidn(t's)  =  1  (5"13) 

and  from  Corollary  5-3, 

Ylup(t's)  ^  y2  u p (t,s)  (5_14) 


for  0  £  s  <_  t. 

The  term  (s)  'Pfj^  (t-s)  =0}  represents  the  probability 
density  function  that  there  are  exactly  l  complete  cycles  in  [0,t] 
and  the  X.th  terminal  point  of  a  cycle  occurs  at  s.  Thus  in 
the  interval  [s,t]  there  is  no  complete  cycle.  Given  this,  the 
probability  that  the  system  will  be  found  in  state  "up"  at 
time  t  is  F^(t-s)  and  the  probability  that  the  process  will 
be  found  in  state  "down"  (under  replacement)  is  [f^G]  (t-s)  . 
Continuing  the  proof  of  Corollary  5-4,  we  rewrite  Eq.  (5-12) 
as 


t  n  ... 

P(N.  (t)  <  n)  -  /  [  l  h,w  (s)  -P(N.  (t-s)  -  0 >  [y .  (t,s) 

1  ~  s-0  k*0  1  1  1  up 

+  y.  dn(t,s))ds 
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-•  >  > 


(5-15) 


p{Ni(t)  <_n} 


A  I  hfk)  (s) ] *P(N. (t-s)  =  0}ds 
0  k=0  1  x 


Then,  from  (5-10)  and  (5-15)  we  get 
t  n  ... 

/  [  l  h-'K'  (s)  ]  •P(N1  (t-s)  =  0}ds 
0  k=0  x  x 

t  n  / 1_  \ 

>  /  [  I  (s)] -P{N.(t-s)  = 0}ds  . 

0  k=0  £  * 

Multiplying  both  sides  of  this  relationship  by  a  non-negative 
function  does  not  change  the  inequality,  so  we  find  that 

=>  A  l  h^k)  (s)]  *p{N_(t-s)  =»  0}-y_  (t,s)ds 

0  k=0  1  1  ^up 

t  n  ... 

<  j  [  l  h'  1  (s) ] *P{N.  (t-s)  « 0}  (t,s)ds 

~  o  k*0  1  x  z  up 

t  n  / \ 

<  J  [  I  h,'K'  (s)  ]  *P(N.  (t-s)  =  0}  ‘Y-.  ,  (t,s)ds 

0  k=0  x  x  x  u  p 

=  A^n) (t)  . 

Q.E.D. 

Corollary  5-5. 

Let  A^,m)  (t)  (i  *  1,  4,  5,  p  *  1,2)  be  the  availa- 
P 

bility  at  time  t  of  a  two-component  system  under  scenario  i 
when  the  system  starts  to  operate  from  state  Dp  with  (n,m) 
spares  available  at  time  0.  Then 

A0i,m)  (t)  < 

A^n#nt)  (t)  « 


A^n,m)  (t) 


A<£'m)  (t) 


(n,m) 


(t)  A 


P 

(n,m) 

1D_ 


<  a(n'm) 
"  1dp 


and 


A^,m)  (t),  p 


1,2 


A*5J  (t) 


•  >  > 
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Proof  of  Corollary  5-5. 
For  the  proof  of 


AU,m)  (t)  <  (t) 


A^n,m)  (t)  <  A^'m)  (t) 


we  let 


fi(t)  "  Vi« 


-B.X.t 


under  scenario  5  (i.e.f  0.^  >  1,  i  =  1,2).  The  recursive 
equations  for  the  availabilities  are: 


A[n,n0  (t)  =  FjFjU)  +  If  ].F2*AiS'm)  J 


+  ff  p  *AIn»m) i 
lt2*l  A1D2  j (t) 


(t) 


(5-16) 


and 


A^n,m)  (fc)  =  F1F2(t)  +  [f1F2*A^'m)] 


(t) 


+  ff  f  *A  1 

ir2*l  A5D2  1  (t) 


(5-17) 


From  Eq.  (5-16)  and  Eq.  (5-17)  we  have 


* (n,m) »  . (n,m) ...  rf  p  _  « (n,m) i. 

A1  A5  w  l£r  2  VA1Dx  A5Dl  n  (t) 

♦  [f/i*(A‘S;m) -A<g'ml)](t)  (5-18 
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Suppose 


AiS'm> ( *)  -  AiS,m)(*) 


-  5D, 


and 


(n,m) 
'ID. 


then  from  Eq.  (5-18) 


A(n,m)  (t)  Agn'm)  (t)  . 


Conversely  suppose 


A^n,m) (fc)  >  A^n,m) (t) 


t  >  0, 


Then  from  Eq.  (5-18)  at  least  one  of  the  two  terms  of  the  right 
hand  side  must  be  positive.  Thus  it  is  sufficient  to  show  that 


AlS'm) (t)  -  A5D'm) (t)  *=*  AiD'm) (t)  -  A5D'm) (t)  * 


That  is,  if  one  of  the  terms  of  the  right  hand  side  of  Eq.  (5-18) 
is  positive,  then  the  other  term  is  positive  whenever 


A^n ,m) (t)  >  A^n,m) (t)  # 


Assume  without  loss  of  generality  that 


AlD'm) (t)  1  A5D'm) (t)  t  -  0 
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We  know  that  by  conditioning  on  the  first  repair  time  of 
component  2  under  each  scenario. 


a1d'B|  *  "i^Kg2t(£1*g1)  (t) 

+  JoIl52l(£l*9l)(lt)*£l*5l11*AiS2'''”'1>>(t)  (5-l9) 

and 


.  (n,ra) 
^5D_ 


n-1 


l  [{g2t(f1*g1>  (]c)*f1']}*a^ 


(n-k,m-l) 


k-0 


(t) 


n-1 


+  J0ll92[(£1*9l1<k)*fI'5lU’A5D1 


(n-k,m-l) 


(t) 


(5-20) 


From  Eq.  (5-10)  we  have 


"-1  »>  .=  “-1 


J0‘<£i**1,lM*F11(u)  +  kr((f1*g1»<k,*£1*olI(u) 


n’1  °  (k)*=T°  n'1  • 


-  J01(£l*3l>'"'*Fl],u,  +  Jot(£l’9l)0'1”£l*°l1(u)  • 

(5-21) 


O  _ 


Multiplying  both  sides  of  Eq.  (5-21)  by  g2(u)  which  assumes 
positive  value  for  all  values  of  u  >  0  yields 


n_1  <*)*=r  n_1 


L{,2((£i*gi>  l,l,*Fi11(U)  +  A  '*2I<£i*gi><k,*£i*5iJ> 


k=»0 


k-0 


(u) 


i  g«gan£I*»i>  (k>*r°) }  (ul  ♦  "l\g2l  t*i*gx)  0t’*  V5i>  >(u)  • 


(5-22) 
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For  ease  of  writing  let 


(g2!(£l*g1)lk)*f1n(u)  > 

af>  („, 

fg2I(£1*g1)(k)*£i*GiH,u. 

-  af»(u) 

{g2l(f“*gi)<k)*p‘])(u)  - 

b'k> (u) 

t92 1 1*!*9!1 <k>  **l*®l'  * (u) 

■  b<k’<u, 

then  the  expression  (5-22)  can  be  written  as 

Y  a‘k>  (u)  *  Y  a<k>  (u)  >  Y  b<k>  (u)  ♦  Y  »‘k>  (»> 

k-0  1  k=0  z  ~  k=0  1  k=0 


(5-23) 


or 


"il  .'kl  (u>  -  Y  b<k>  (u)  >  Y  b<k>  (u)  -  Y  af»  (u,  . 

k*0  A  k=0  x  k=0  k=0 


(5-24) 


Ey  the  assumption  of 


A(k,£)  (t)  >  (t)  for  all  integers  k,i  >_  0 

1  5  t  >  0. 


and  by  the  inequality  of  (5-24),  convolution  operations 
show 
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n-1 
k 


n-1 


n-1 


^[a^k)*A^n‘k'in‘1)]  (t)  -  kIo[bik)*A^n'k'm"1)  ]  (t) 

i  -  ny1rb(k)*A(n-k,m-: 

J(t)  k=ol  1  1  (t) 

(k) (n-k,m-l) 


(k) *  (n-k,m-l) 


k=0  1  1 


3 


>  ni1tb<k|*A<n-k'”-l)],fl  -  ni1[a‘k)*»1(n-k',"-1)) 

~  k=0  z  1  [,Z)  k=0  z  1 


n-1 


>(k) *, (n-k,m-l) 


l  (*2  *A1D 
k=0  2  iUl 


n-1 


n-1 


(k) ±  (n-k,m-l) 


1  -  7  faliC,*A' 

](t)  k£0[a2  A1Dx 

n-1 


(t) 


(t) 


>  7  fb  ^k^  *A^n-k,m~^  1  -  7  [ a  ^k^  A (n-k*m“3-)  ]  (5- 

-  .i.lD2  A5D.  1  (t)  k£0l  2  ID,  J  (t)  1 


k=o  2  5D1 


(t) 


Rewriting  expression  (5-25)  yields 


n-1 


(k) *. (n-k,m-l) 


l  [a^  *A, 

k=0  1  1 

n-1 


n-1 


(k) *, (n-k,m-l) 


i  +  y  raw*A' 

(t)  k*0  2  1D1 


(t) 


>  l  [b'*'*A* 
~  k=0  1  5 


(k) t  (n-k,m-l) 


n-1 


i  (t) +  i  fbr*A 

k=0  * 


(k) (n-k,m-l) 


5D. 


(t) 


which  we  recognize  as 


, (n,m) ...  « (n,m) 

A1D2  -  A5D2 


from  the  expressions  (5-19)  and  (5-  20>  .  We  have  shown  that 
whenever 


A|n,m)  (t)  >  A^n ,m)  (t) 


holds  .  it  follows  that 


25) 
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AU,nO  (fc)  ^  A^n,m)  (t)  and 

A^(n,m)  ^  A^n,m)  ^  are  eqUivaient 

P  P 

for  p  =  1,2. 

In  a  similar  fashion  it  can  be  shown  that 


A^n,m) (fc)  ^  A^n,m) (t)  and 

A^n,m) ^  ^  A^n,m)  are  eqUivaient 

P  P 

for  p  *  1,2,  which  completes  the  proof  of  Corollary  5-5. 
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Theorem  5-6 


(n,m) 


(n,m) 


(n,m) 


'5  (t)  1  hi  (t)  1  Aj  '  1  (t)  t  >_  0,  m,m  «£  N 


Proof . 

For  the  proof  of 


A^n,m) (t)  ^  A|n,m) (t) 


Let 


O 

fi(t) 


"9iXit 

9iXie  A  for  0i  >  1  ,  i  ■  1,2 , 


and  i  =  1,2.  When  (n,m)  *  (0,0) 


A<°'0,<t>  =  A<°'0)(t)  - 


-Uj+Xjlt 


and  when  (n,m)  =  (1,0) 


A;1-°»(t) 


FlF2(t>  +  (t)  i  A^1'0*  (t) 


F1F2  +  [£lt'2*9lF2*FlF21  (t) 


The  last  Inequality  holds  because 


[9lF^(t)  i  [3lF2I 


(t) 


T  >  0  . 
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Liz. 


Also,  when  (n,m)  =  (0,1)  we  have 


A^0,1)  (t)  <_  A^0,1)  (t) 


We  will  assume  that 


A^1,j)(t)  1  A^1,j)(t) 


for  (i,j)  <  (n,m)  and  show  that 


AU,m)  (t)  ^  A^n,m)  (fc) 


From  Eq .  (5-5)  we  have 


A(n,m)(t)  _  PlP2(t)  +  ”j\f1Ef(g1[(£2*g2)  (k,*F2l)l 


A<n-l-n.-k)(t-sl)dsl 


t  m-1 


+  S0  ltI0l£l^<glI<£2*g2)<k)*,2*G21H(  =  1) 


*  J/£2ri*<92l(£l*5l»“)*Fll>l(s,) 


A|n-i,m-l) (t_S2)dS2 


♦  ,,t”yf^i,2|(ti*,i,l,l’'i‘0iin'*2' 


tv  ( n* &  f in- 1)  /  l  \  j_ 

•  AlDl  (t-s2)ds2. 


In  a  similar  fashion,  we  set  up 
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A^n,m) (fc) 


a>k)(3)  =  (f1P2*(gi[(£2*g2)(k>*F2l,](s) 

b<k)(sl  =  [£lV<9lI(£2*92,(k’*£2*521>)(s) 

a‘e,(s)  =  t£2F1*Ig2[(£1*g1)(t)*F1l))(s) 

b2t’<s>  *  t£2Fl*I52[I£l*9l|(‘,*£l*5lI,1(s) 
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and  we  let  a^  f  b^k^  (s) ,  (s)  /  b^  (s)  be  the  corresponding 

terms  in  Eq.  (5-27). 

To  prove  that 

A^n'm)  (t)  _>  A^n,m)  (t) 

we  must  show  that 

my1{ ra(k)*A(n_1,m“k) 1  +  [b(k)*A(n-1'm-k)l  } 

k=0  1  1  (t)  1  1D2  J (t) 

>  my1{ra°(k)*A(n"1'm‘k)]  +  rb°(k)*A(n-l,m-k) 

-  ,1^  lal  A5  J  (t)  +  li3l  A5D„  1  (t)  } 


k=Q 


(5-28) 


and 


n-1 


y  {  Ta  ^  1  +  fb  ^ ^  *A ^n-Jl  1  } 

L  la2  A1  ]  (t)  +  it32  A1D.  J  (t)  * 


1=0 


n-1 


i  J.{  la2(U*A5n~li'ro"1)  J  ( t)  +  m  > 


1=0 


5D, 


(t) 
(5-29) 


This  is  accomplished  beginning  with  expression  (5-10)  which 
implies  that 


B£j<V9,)(k)*F2l(u)  +  ”iVf2*g2)(lt)*£,*<=2] 

U 


k=0 


i  J  niV92l<k,<l,..,  +  t(£,*g,)',1'*£;*G,),,il} 


k=0 


2 J  (u) 


O 

'2  *2 


2  2  (u) 
00**°*^ 


■2  2  (p) 


(5-3  0) 
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1 


Multiplying  this  expression  by  g^(u)  and  taking  convolution 
by  f^F^(-)  on  both  sides  of  (5-30)  yields 


m-1 


*0  1  lsl’  k-0  1  l8l'  k-0  1  lsl' 


m-1 


(k) 


m-1 


(k) 


,  my*  iw°  (k) , 

Jo'1  <S1> 


(5-31) 


Since 


A|i,j)(*)  >  A^'^t*) 


for  (i,j)  <  (n,m) ,  convolution  operations  show  that 


m-1 

l 

k=0 


y  /  r a  (k)  * j.  (n—  1  /Hi— k)  i  _  r a°  (k)  *»  (n—  1  »m—  k) .  . 

£ftUal  A1  J (t)  lal  A5  J(t)} 

y  r  r-.(k)  (n-l,m-k) ,  _  ,  °  (k)  ^  (n-lf  m-k) ,  , 

kio  1  1  1  (t)  1  1  J  (t) 

y  r  rjj°  (k)  * .  (n— 1 ,m— k) ,  ,,  (k)  #^(n-l,m-k) ,  , 

-  k£0l lbi  Ai  J (t)  lbi  Ai  J (t) * 

y  /  [Vj°  ^k)  (n-l,m-k) ,  _  ,^(k)*  (n-l/m-k)j  . 

-  k£0  i  aid2  J(t)  lDi  aid2  J  (tr 

y  /  rvj°  (k)  (n-l,m-k)  ,  [h  (k)  *&  (n-1  ,m-k) ,  , 

-  l“l  A5D2  J  lDl  A1D2  J  (t)  1  ‘ 

The  second  inequality  above  comes  from  (5-31),  the  third 
inequality  comes  from 


*{i,31<->  i  a^;3,(-> 
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and  the  last  inequality  comes  from  Corollary  (5-5)  .  Thus, 
we  have  verified  (5- 28) .  In  a  similar  fashion  we  cam  verify 
the  expression  (5-29).  Thus,  we  have  proved 

A(n,m) ^  A^n ,m) (t) 
for  t  _>  0  and  n,m  €  K. 

Following  the  exact  procedures  used  in  the  proof  of 
A^n,m) (t)  ^  A^n,m) (t) 

it  can  be  shown  that  the  inequality 

Ajn,m) (t)  ^  A^n'm) (t) 

is  true,  proving  Theorem  5-6. 

Corollary  5-7 : 


(n,m) 

5dp 


it) 


< 


a  (n,m) 

1dp 


it) 


a (n,m) 

4dp 


(t>  . 


From  Theorem  5-6  we  have  Agn,m^  (t)  <_  A|n,m^  (t)  .  By  Corollaury 
5-5  this  implies  that  <_  A^'m^  (t)  .  Also,  we  have 

A^n.m)  (t)  ^  A^n,m)  (t)  which  implies^" ,m)  (t)  <_  A^'m)  (t)  . 


a (n,m) 
a5D 

P 


(t) 


< 


* (n,m) 

10p 


it) 


a (n,m) 

4Dp 
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VI.  ALLOCATION  OF  SPARE  PARTS 


A.  BACKGROUND 

In  the  previous  chapters  we  have  defined  operational  availa¬ 
bility;  we  have  derived  computational  formulae  for  the  point 
availability  when  only  finitely  many  spare  parts  are  available; 
we  have  looked  at  the  impact  of  the  operational  deployment 
scenario  and  configuration  on  system  point  availability ;  and 
we  have  determined  order  relationships  for  system  availability 
under  different  scenarios.  The  original  motivation  for  this 
research  effort  was  to  develop  a  methodology  for  allocating 
scarce  resources  (in  this  problem,  dollars  spent  for  spares) 
so  as  to  maximize  system  availability.  In  this  chapter  we 
develop  such  an  algorithm  which  utilizes  the  results  obtained 
in  the  earlier  chapters  to  optimally  allocate  dollars  for 
spare  parts  for  the  components  within  a  given  system.  Depending 
on  one's  definition  of  the  "systan",  this  algorithm  can  be  applied 
toward  the  allocation  of  spare  parts  for,  say,  an  entire  ship 
or  for  the  individual  components  of  a  weapon  such  as  a  radar 
or  a  fire  control  system. 

Various  authors  have  considered  resource  allocation  of 
spares  (parallel  redundancy)  in  order  to  maximize  system  relia¬ 
bility  with  a  budget  constraint.  Moskowitz  and  McLean  [30] 
used  a  variational  method  to  optimize  redundancy.  Burton  and 
Howard  [  7 ]  and  Bellman  and  Dreyfus  [  6 ]  solved  the  problem  using 
dynamic  programming.  Kettelle  [22]  developed  a  heuristic  algorithm 
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from  dynamic  programming  and  that  was  extended  by  Proschan 
and  Bray  [33]  to  allow  for  multiple  constraints.  c’lare  and 
Taylor  [18]  used  a  branch-and-bound  procedure;  Mizukami  [29] 
used  convex  integer  programming,  and  Tillman  and  Liittschwager  [44] 
maximized  reliability  subject  to  several  constraints  using  am 
integer  programming  formulation.  Everett  [14]  used  the  generalized 
Lagrange  multiplier  method  to  solve  the  same  problem  discussed 
by  Kettelle.  Sharma  and  Venkateswaran  [38]  and  Nakagawa  and 
Nakashima  [31]  developed  intuitive  algorithms  that  provide 
approximate  solutions .  In  all  of  those  cases  the  system  con¬ 
sists  of  components  in  series  with  parallel  redundancy.  Aggarwal 
[  1  ] ,  and  Kuo,  Hwang,  and  Tillman  [24]  develop  heuristic  methods 
for  optimal  system  reliability  for  more  complex  systems.  Other 
papers  [39]  and  [40]  have  attempted  to  determine  stock  levels 
for  components  in  order  to  maximize  equipment  operational 
availability  ,  subject  to  a  budget  constraint.  For  their 
purposes  operational  availability  for  an  equipment  is  calculated 
using  the  definition 

.  _ MTBF _ 

A  MTBF  +  MTTR  +  MSRT 

where  MSRT  is  the  mean  supply  response  time.  Since  MTBF  and 
MTTR  are  unaffected  by  the  number  of  spares,  the  models  really 
minimize  mean  supply  response  time.  The  method  they  use  is  a 
Lagrange  multiplier  approach  with  an  embedded  dynamic  programming 
technique . 
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In  the  research  summarized  above  there  have  been  many 
approaches  to  solving  resource  allocation  problems  using  vari¬ 
ous  mathematical  programming  techniques  and  considering  differ¬ 


ent  types  of  resource  constraints.  Our  spares  allocation 
problem  could  also  be  solved  several  different  ways.  The 
system  availability  objective  function  does  not  possess  the 
required  characteristics  (linear,  convex,  concave)  that  allow 
the  use  of  some  of  the  specialized  programming  techniques. 

It  is,  however,  separable  and  monotone  in  the  decision  varia¬ 
bles  (the  numbers  of  spares)  and  the  constraints  are  linear. 
Therefore,  solution  techniques  like  dynamic  programming  or 
heuristic  methods  can  be  applied. 

None  of  the  algorithms  described  above  were  developed  for 
the  availability  allocation  problem  that  we  address.  However, 
there  are  similarities  that  we  can  exploit.  In  the  next 
section  we  present  an  allocation  algorithm  for  the  availability 
problem  that  extends  the  results  obtained  by  Burton  and  Howard 
[7]  and  Kettelle  [22] . 

B.  THE  ALLOCATION  ALGORITHM 

Consider  a  weapon  system  with  k  components,  and  let  n^  be 
the  number  of  spares  allocated  to  component  i  for  i  =  l,2,...,k. 
Let  L^  and  be  lower  and  upper  bounds  on  the  number  of  spares 
for  component  i  and  let  n  =  (n^,n2, . • . fn^)  be  the  vector  of 
spares  allocation  for  the  k  components.  Finally,  let  c (n)  be 
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the  total  cost  for  the  allocation  n  and  B  the  upper  bound  on 
the  dollars  available  for  allocation.  The  mathematical  prob¬ 
lem  we  address  is 


max  A^  (t) 
subject  to  c  (n)  <_  B 


i  ■  1,2, ...  ,k 


(6-1) 


where  A^— ^  (t)  is  the  system  point  availability  with  allocation 
n. 

The  mathematical  programming  problem  (6-1)  has  little  in 
the  way  of  special  structure  in  the  objective  function  which 
allows  the  use  of  various  efficient  linear  or  nonlinear  pro¬ 
gramming  algorithms.  However  the  system  availability  does 
possess  a  separability  and  monotonic  structure  that  allows 
the  use  of  dynamic  programming  methods.  The  solution  technique 
that  we  derive  is  based  on  Kettelle's  reliability  redundancy 
allocation  algorithm  [22]. 

Kettelle  provides  an  easily  usable  algorithm  for  obtaining 
an  exact  solution  for  maximizing  the  reliability  of  a  parallel 
redundancy  series  system  subject  to  a  budget  constraint.  His 
algorithm  generates  undominated  redundancy  allocations  (or 
dominating  sequence  of  allocations)  for  successively  larger 
subsystems  from  undominated  allocations  for  small  subsystems. 


To  understand  the  concept  of  "an  undominated  allocation" 


the  following  definition  is  introduced  [5] . 

Definition:  n°  is  undominated  (or  dominating)  if 

A ^  (t)  >  A^—  ^ (t)  implies  c(n)  >  c(n°) 


whereas 


A^n^ (t)  =  A^—  ^ (t)  implies  either  c(n)  >  c(n°) 


where 


or  c(n)  =  c(n°) 


Kettelle  confined  his  algorithm  to  series-type  systems. 

Burton  &  Howard  [  7 ]  showed  that  the  allocation  problem  can 
be  solved  by  dynamic  programming  for  any  system  configuration 
composed  of  a  mixture  of  series  and  parallel  connections.  They 
developed  a  computer  algorithm  for  this  problem  and  demonstrated 
that  the  dynamic  programming  method  works  very  well  for  complex 
sys  terns . 

The  Burton-Howard  recursive  dynamic  programming  algorithm 
can  be  improved  computationally  by  introducing  Kettelle* s  idea 
of  undominated  sequences  and  his  tableau  computation  methods. 

The  algorithm  which  we  develop  for  optimizing  system  availability 
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adapts  features  of  both  the  Kettelle  algorithm  and  the  Burton- 
Howard  algorithm  and  results  in  an  algorithm  which  is  compu¬ 
tationally  efficient.  The  result  solves  the  allocation  problem 
not  only  for  the  budget  B,  but  also  for  all  budgets  B1  <_  B. 

This  is  important  since  in  real  world  applications  the  budget 
B  itself  is  not  always  known  precisely  at  the  time  of  solution. 

We  describe  the  algorithm  in  the  following  material.  We 
assume  that  the  computational  formulae  derived  earlier  are  used 
to  evaluate  component  availability  and  that  the  system  con¬ 
figuration  is  known.  A  numerical  example  is  provided  in  Section 
VI. C. 

The  steps  of  the  availability  allocation  algorithm  are 
described  below.  Since  the  algorithm  was  generated  out  of  a 
dynamic  programming  solution,  it  is  not  necessary  to  prove 
optimality.  (Proofs  are  available  from  the  Kettelle  and  Burton- 
Howard  references.)  Where  we  have  made  modifications  of  the 
previous  results  to  reduce  the  number  of  decision  alternatives 
at  a  given  point,  we  prove  that  the  modifications  cannot  result 
in  inferior  solutions. 

The  Algorithm 

Define  the  N  stages  (N  <_  k)  to  consist  of  the  independent 

system  entities.  These  entities  can  be  the  individual  components 

if  they  operate  independently  or  modules  composed  of  dependently 

operating  components  (such  as  two  components  in  series  operating 

(n  > 

under  scenario  2) .  Let  a^  be  the  availability  (for  con¬ 
venience  we  drop  reference  to  time  t  even  though  there  is  still 
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a  specific  time  t  in  which  we  are  interested)  for  the  indepen¬ 
dent  entity  i,  i  =  1,2,...,N. 

(n  ) 

1)  Compute  the  stage  return,  a^  for  stage  1  for 

Li  <_  n^  <_  U^.  This  calculation  utilizes  the  results  derived 
in  Chapters  II  through  V. 

2)  Formulate  the  N-stage  return  function  using  the  availa¬ 
bility  calculus  for  independent  series/parallel  configurations. 
Burton  and  Howard  show  that  the  problem  can  always  be  formulated 
so  that  separability  and  monotonicity  are  satisfied  for  any 
series/parallel  mixed  system  configuration  provided  the  entities 
are  independent. 

3)  Set  up  a  tableau  such  as  that  shown  in  Tables  6.7 

for  the  two-stage  problem.  The  entries  in  the  row  headings 

(  1> 

correspond  to  the  triple,  (n^n^c^,  and  a^  );  the  entries 

( n  _ ) 

in  the  column  heading  correspond  to  (r>2,n2c2,a2  ).  The 

entries  in  the  body  of  the  table  give  the  allocation  (n1,n2>, 
the  cost  cfn^n^  and  the  two-stage  return  (availability). 

4)  Start  with  (1,1)  as  the  first  element  in  the  sequence 
of  undominated  allocations  for  the  subsystem  consisting  of 
s  tages  one  and  two . 

5)  Select  as  the  next  member  of  the  sequence  of  dominating 
allocations  the  cheapest  cost  entry  with  availability  higher 
than  the  previous  element  of  the  sequence.  (We  discuss  later 
some  methods  of  eliminating  entire  blocks  of  possibilities 
from  consideration.) 

6)  After  proceeding  through  a  given  tableau  in  the  above 
manner,  increase  the  problem  to  a  3-stage  return  using  the 
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undominated  allocations  from  the  previous  2-stage  problem 

(n.) 

as  the  row  entries  and  (n3,c3 (n3) ,a3  J  )  as  the  column  entries. 

Obtain  the  dominating  allocations  for  this  subsystem  consisting 

of  stages  1,  2,  and  3  as  before. 

7)  In  general,  for  the  N-stage  problem,  use  as  the  row 

(n, , n_ , . . • , n„. ) 

elements  ( (n^,n^, . • • , ,c(n^,n2,..», ^ )  ,  A  ) 

and  as  the  column  elements  <nN,c(nN),a  )  and  step  through 
selecting  a  sequence  of  dominating  allocations  by  taking  the 
cheapest  cost  entry  with  availability  higher  than  the  previous 
value . 

Note  that  if  c^  and  a^  represent  the  cost  and  availability 
of  cell  (i,j),  and  if  cell  (i,j)  is  a  member  of  the  optimal 
sequence,  then  the  next  member  of  the  optimal  sequence  must 
be  located  in  the  region  represented  as  quadrant  I  or  III  in 
Figure  6.1  below. 


Fig.  6.1  Definition  of  quadrant  I,  II,  III,  IV 

The  tableau  method  used  in  this  algorithm  has  the  following 
advantages : 

(1)  It  does  not  need  backtracking  procedures  which  are  used 


in  the  standard  dynamic  programming  algorithm. 


(2)  It  does  not  need  extra  calculations  for  the  changes  of 
budget  level.  Since  we  can  read  off  Figure  6.3,  the  allocation 
and  availability  for  each  level  of  the  budget  ranges  from  0 

to  53. 

(3)  It  is  efficient  especially  for  manual  calculations.  As 
we  could  see  in  the  previous  tables  (Tables  6.7  -6.10)  we 

do  not  need  to  fill  out  every  cell  of  the  tables.  Heuristic 
inspection  of  some  of  the  cells  allows  the  user  to  eliminate 
large  blocks  of  cells.  As  an  example,  suppose,  in  Table  6.9, 
we  have  just  computed  elements  of  the  cell  (9,7)  of  the  table 
6.9  which  shows  allocation  (3, 4, 0,0, 8),  cost  17.4  and  availa¬ 
bility  0.6336,  and  we  found  out  that  this  allocation  should 
currently  be  included  in  the  undominated  sequence  according  to 
the  5th  criterion  of  this  algorithm,  then  we  know  that  we  have 
already  considered  all  the  elements  of  the  cells  above  the 
solid  line  in  Table  6.9,  all  of  which  have  costs  less  than  17.4. 
So  in  the  next  calculations,  we  may  try  the  entries  in  several 
cells  in  quadrant  I  or  III  for  the  possible  elimination  of 
entire  blocks.  Suppose  for  example  we  compute  the  availability 
in  the  ceil  (19,5)  of  Table  6.9  which  we  found  to  be  0.5966 
(lower  than  0.6336)  with  cost  24.6  (higher  than  17.4),  then 
we  do  not  need  to  consider  cells  which  are  contained  in  the 
shaded  area  ("North-West  Corner"  elimination  rule) .  This  is 
because  cost  and  availability  of  the  cell  (19,5)  of  Table  6.9 
are  maximal  among  those  of  the  cells  contained  in  the  shaded 
area  (cost  and  availability  in  the  cell  (i,j)  of  the  table  are 


134 


increasing  in  i  and/or  j),  and  thus  every  cell  in  the  shaded 
area  has  lower  availability  with  higher  cost  than  whose  of 
cell  (9,7)  which  is  the  current  element  in  the  dominating 
sequence . 

In  those  stages  where  cumulative  stage  returns  are  com¬ 
puted  from  series  connections,  the  remainder  of  an  entire  row 
or  column  can  sometimes  be  rejected  as  dominated.  Let  c^q 
and  a^Q  represent  the  cost  and  availability  of  the  entry 
heading  row  i.  If  >  a^fQ  where  i'  <  i,  all  entires  in 

row  i'  which  cost  more  than  c^j  are  dominated  because 
a. .  >  a. ,  >  a. ,,  k  >  1.  The  second  inequality  follows  because 

IO  1  K  “■ 

a^,k  =  a^i0*ac>jc  where  aQk  <  1,  k  =  1,2,...  .  The  same  is 
true  for  columns.  As  an  example,  note  that  in  Table  6.7 
cell  (2,2)  (cost  5.4,  availability  0.330135)  dominates  every 
cell  in  the  first  row  starting  with  (1,4)  since  a22  =  0.330135 
>  alo  =  0.2865  >  a^^  (k  >_  1)  ,  and  c22  =  5.4  <  c^£,  l  >_  3 . 

In  those  stages  where  cumulative  stage  returns  are  com¬ 
puted  from  parallel  connections  there  seems  to  be  no  clear 
cut  rule  of  elimination  other  than  the  "north-west  corner 
elimination  rule",  i.e.,  try  several  cells  heuristically  in 
quadrant  I  and  III,  if  a^,j,  <  then  we  eliminate  from 

further  consideration  cell  (n,m)  such  that  1  <_  i'  and  m  <_  j  ' . 

C .  AN  EXAMPLE 

As  an  illustration  consider  the  system  configuration  shown 
in  Figure  6.2  and  the  data  in  Table  6.1.  Solve  the  following 
problem. 


1 


The  problem  is: 


max  A^  <t) 


s.t.  n  •  c  <  B 


0  <  n.  <  U. 


for  t  =  100,  B  is  flexible  between  40  and  50. 

Solution:  For  this  problem  the  availability  is  given  by 


(nl'n2'n3 ,n4 ,n56'n78* 


(n7Q)  (n.^.) 

1-  [l-a?878  ][l-a5656  ] 


(n7,n  )  (n  )  (n_) 

{1  -  (1  -  a34J  q  )  (1  -  a  1  a2  z  ) }] 


(6-2) 


which  is  separable  since  it  can  be  written  as 


(n. »no »n  ,n . ,n_fi ,n7R ) 

A  "  a78  °  a56  °  a34  °  a2  °  ai 


where  "o"  represents  the  composition  operator. 
We  seek 


=  max  A 


(ni'n2 ,n3 ,n4 ,n56 ,n78^ 


s.t.  n  •  c  <  B 


0  !  n.  <  U.  . 


Let  a.^  be  the  stage  return  and  A^(X^)  be  the  maximum  i-stage 


return. 
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For  stage  1, 


(iO 

A1*X1^  =  max  al 

nl 


s.t.  0  <  nj  <  6 

(6-3) 

11  1  cini  1 

For  stages  2,  3,  4 

and  5 

a2(x 

(n-) 

2)  =  max  a2  .A.  (X2~n2c2) 

n2 

s.t.  0  <  n2  <  6 

(6-4) 

0  <_  n2c2  <  X2 

A34(X34}  =  max 

n3'n4 

(n3,n.) 

(1  -  [1  -  a34J  4  ] [1  -  A2(X34  - 

n,c_  -  n.c 
3  3  4 

s .  t. 

0  i  n3  <  3 

0  i  n4  1  n5 

(6-5) 

n3C3  +  n4C4  -  X34 

AecfX-,)  =  max 

56  56  "56 

(n565 

a56  ,A34  (X56  '  C56n56) 

s .  t. 

0  i  n56  ^  6 

(6-6) 

0  1  c56n56  i  X56 
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*78^*78*  =  1  “  min  ^1-a78  "  A56  ^X78  **  C78n78^ 

n7  8 

s-t-  0ln78l7  (6-7) 

0  -  C78n78  -  X78  =  B 
(n  ) 

The  stage  returns  (a^  )  are  computed  from  the  formulas  in 

Chapter  IV  and  the  computation  results  are  listed  in  Tables 
(6.2  -  6.6) .  Here  n^  is  the  total  number  of  parts  (original 
plus  spares  for  component  i) . 

At  stage  1  and  2,  from  Eq.  (6-3)  and  Table  6.2  the  maximum 
return  from  stage  1  is 


A1(n1c1) 


(n^) 

al 


To  obtain  a  complete  sequence  of  undominated  allocations  for 
the  subsystem  consisting  of  stages  1  and  2  according  to 
recursive  Equation  (6-4),  we  set  up  Table  6.7.  The  entires 
in  the  body  of  the  Table  (6.5)  give  the  spares,  (vector  of 
two  elements) ,  cost,  and  availability  for  the  subsystem  con¬ 
sisting  of  stages  1  and  2.  Thus,  the  entry  (2,3)  corresponds 
to  n^  =  2  and  n£  =  3  with  cost  6.7  achieving  subsystem  availa¬ 
bility  0.495725.  The  chosen  elements  connected  by  arrows  form 
an  undominated  sequence  of  allocations . 

The  elements  of  the  sequence  are  chosen  in  the  following 

way. 

(1)  Start  with  (1.1),  the  first  undominated  allocation. 

(2)  The  next  undominated  allocation  is  the  cheapest  cost 
entry  with  availability  higher  than  that  of  previous  allocations 
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These  procedures  are  consistent  with  the  Burton-Howard 
algorithm.  The  sequence  of  undominated  allocation  in  Table 
6.7  shows  that  if  we  have  15  as  our  budget  for  this  two- 
component  subsystem  then  we  allocate  n^  *  5  and  n2  =  6  with 
cost  14.8  achieving  availability  0.89055,  if  B  *  20  then 
n  =  (7,7)  with  cost  19.9  achieving  0.8964  and  so  forth. 

At  stage  3,  component  #3  and  #4  form  a  two-component  series 
subsystem  which  operates  under  scenario  2  with  spares  not 
shared.  The  stage  returns  are  computed  from  the  formula  (4-22) 
and  (4-23)  and  listed  in  Table  6.4.  The  sequence  of  undominated 
allocations  is  marked  by  arrows  also  in  Table  6.4.  From 
recursive  equation  (6-5)  we  obtain  a  sequence  of  undominated 
allocations  for  the  subsystem  consisting  of  component  #1,  #2, 

#3  and  #4  by  using  the  undominated  allocations  of  stage  1,  2 
and  the  undominated  allocations  of  stage  3  which  is  shown  in 
Table  6.8. 


At  stage  4,  since  component  #5  and  #6  operate  under  scenario 

2  with  spares  shared  by  both  components,  the  payoff  function 
(n56> 

a56  whihc  is  computed  from  formula  (4-27)  and  (4-28)  and 
listed  in  Table  6.5  is  the  availability  of  the  module  which 
is  an  independent  entity. 

According  to  the  recursive  equations  (6-6)  we  obtain  a 
sequence  of  undominated  allocations  for  the  subsystem  con¬ 
sisting  of  components  #1,  #2,  #3,  #4,  #5,  and  #6  by  using  the 
undominated  allocations  of  stages  1,  2,  3  and  the  sequence  of 
payoffs  from  stage  4  (component  #5,  #6).  Table  6.9  shows  this 
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sequence  ,  and  Table  6.10  shows  final  allocations  according 
to  the  recursive  equations  (6-7)  in  a  similar  fashion. 

From  the  final  table  (6.10)  we  can  construct  the  following 
graph  which  shows  allocation  and  availability  for  each  level 
of  budget  which  is  flexible  ranging  from  40  to  50. 

If  we  would  like  to  have  a  system  availability  of  0.900 
then  we  would  need  only  a  budget  of  26.7  consumed  by  allocating 
(n1,n2,n3,n4,n56,n7g)  =  (3 , 3 , 0,0 , 6 , 9) .  Figure  6.3  is  the 
graphical  representation  of  Table  6.10. 


Fig.  6.3  Availability-cost  curve  for  undominated 
allocations  (budget  level  40  -  50) . 
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Table  6-10:  Sequence  of  Max  Return  from  Stage  1,  2f  3,  4, 
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VII.  SUMMARY  AND  CONCLUSIONS 

One  of  the  biggest  problems  facing  the  military  services 
is  that  of  measuring  the  readiness  condition  of  their  fighting 
units.  In  this  thesis  we  have  been  concerned  only  with  the 
material  readiness  of  weapon  systems.  The  material  readiness 
of  a  weapon  system  depends  on  many  factors  among  which  are 
equipment  reliability,  spare  parts  support,  supply  response 
time,  maintenance  philosophy,  maintenance  manning  levels,  sys¬ 
tem  configuration,  and  operational  deployment  scenario.  Many 
parochial  measures  of  material  readiness  are  calculated  and 
reported  by  separate  organizations  within  the  services.  For 
example,  the  supply  organizations  calculate  indicators  of  the 
level  of  supply  support  such  as  fill  rate  or  supply  response 
times  while  maintenance  organizations  report  on  repair  times, 
and  others  are  interested  in  reliability.  One  widely  reported 
measure  of  material  readiness  which  attempts  to  take  into 
account  all  of  the  avoe  listed  factors  is  operational  availa¬ 
bility.  Because  operational  availability  considers  a  multi¬ 
tude  of  factors,  it  is  the  type  of  measure  on  which  higher 
command  levels  tend  to  focus. 

Operational  availability  is  widely  understood  to  be  a 
measure  of  the  likelihood  that  a  system  will  function  success¬ 
fully  when  called  upon.  However,  there  are  differences  in  the 
way  that  operational  availability  is  calculated.  Military 
documents  specify  that  operational  availability  be  calculated 
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by  taking  the  ratio  between  the  mean  time  to  failure  and  the 
sum  of  the  mean  time  to  failure  and  the  mean  time  to  replace 
a  failed  unit.  This  definition  ignores  many  of  the  critical 
factors  of  interest  and  it  is  often  not  mathematically  correct. 
Other  definitions  of  operational  availability  focus  on  the 
probability  that  the  system  will  be  operational  at  a  specified 
point  in  time  (the  point  availability) ,  the  integral  of  point 
availability  over  a  given  period  of  time  (the  average  availa¬ 
bility)  ,  and  the  limiting  probability  that  the  system  will  be 
available  (steady  state  availability) .  In  this  report  we  have 
concentrated  on  the  point  availability  since  each  of  the  other 
two  availability  expressions  can  be  found  from  point  availability. 

The  primary  objective  of  this  research  was  to  develop  a 
procedure  which  will  optimally  allocate  spare  parts  to  the 
components  of  a  weapon  system  so  as  to  maximize  the  point  availa¬ 
bility  of  the  system  within  a  budget  limit  for  any  specified 
point  in  time.  The  final  result  of  this  thesis  is  an  algorithm 
which  accomplishes  that  objective. 

The  weapon  system  implied  in  this  study  is  composed  of 
components  connected  in  series  and  parallel.  Each  component 
is  subject  to  random  failures  and  random  restoration  times  until 
all  replacement  spares  are  exhausted.  The  failure  times  of 
each  component  are  assumed  to  be  independent  and  identically 
distributed  with  exponential  distributions.  Similarly,  the 
restoration  times  for  a  component  are  exponential  IID  random 
variables.  We  have  considered  cases  where  the  components 
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operate  independently  of  each  other  and  cases  where  the  com¬ 
ponents  do  not  operate  independently.  For  those  systems  with 
components  which  all  operate  independently,  the  system  availa¬ 
bility  can  be  determined  easily  after  the  point  availability 
of  each  component  has  been  determined  as  a  function  of  the 
number  of  spares  available. 

For  systems  containing  components  which  do  not  operate 
independently,  the  system  point  availability  cannot  be  deter¬ 
mined  from  an  application  of  the  standard  reliability  calculus 
of  series/parallel/mixed  systems.  Instead,  special  considera¬ 
tion  must  be  given  to  the  type  of  dependence  that  exists  among 
the  components.  Five  different  operating  scenarios  are  con¬ 
sidered  for  a  two  component  series  system.  These  operating 
scenarios  are  distinguished  by  what  happens  to  the  survivors 
while  a  failed  component  is  being  replaced.  Under  scenario  one, 
the  non-f ailed  component  continues  in  service,  and  if  it  fails, 
its  replacement  begins  immediately  and  proceeds  independently 
(the  independent  scenario) .  Under  scenario  two,  the  surviving 
component  is  shut  down  until  replacement  of  the  failed  component 
is  accomplished  (the  symmetric  shut-down  scenario) .  Under 
scenario  three,  a  failure  of  component  number  one  does  not 
affect  the  operation  of  component  number  two,  but  a  failure 
of  component  two  will  shut  down  component  number  one  (the 
asymmetric  shut-down  scenario) .  For  scenarios  four  and  five, 
the  failure  of  one  component  causes  the  failure  rate  of  the 
survivor  to  change  from  A,  to  A.0..  In  scenario  four  0.  <  1,  and 
in  scenario  five  0^  >  1. 
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Exact  expressions  for  the  point  availability  of  a  one 
component  system  and  for  the  marginal  contribution  to  availa¬ 
bility  of  the  nth  spare  are  determined  as  a  function  of  the 
failure  and  repair  parameters  and  the  number  of  spares.  For 
large  values  of  n,  normal  approximation  are  obtained  to  provide 
simpler  expressions. 

The  exact  expression  for  the  average  availability  is  deter¬ 
mined  from  the  point  availability.  The  tradeoff  between  pro¬ 
viding  a  maintenance  capability  (implicitly  equivalent  to  an 
assumption  of  infinitely  many  spares  and  larger  replacement 
times)  and  providing  modular  replacement  with  finitely  many 
spares  is  discussed .  The  maintenance  repair  rate  n '  that 
achieves  the  same  availability  provided  by  the  replacement 
policy  with  n  spares  is  determined.  When  n  is  small,  the  re¬ 
placement  rate  n  must  be  very  large  compared  to  the  repair 
rate  n ' ;  however,  when  n  gets  large  the  ratio  of  n '  to  n 
approaches  one. 

Given  a  fixed  length  of  mission  duration  and  a  finite 
number  of  spares  a  system  may  not  be  available  at  the  end  of 
a  mission  due  to  a  lack  of  spares.  The  probability  distribution 
of  this  downtime  is  determined  as  a  function  of  the  number  of 
spares  and  mission  duration.  This  probability  distribution  could 
be  used  to  generate  other  measures  of  effectiveness  for  deter¬ 
mining  the  number  of  spares  that  should  be  allocated  to  a  system. 

For  systems  containing  two  or  more  components,  the  system 
availability  can  be  determined  easily  by  calculations  identical 
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to  those  used  to  calculate  system  reliability  from  component 
reliabilities  whenever  the  components  operate  independently. 

All  that  is  needed  is  to  replace  component  reliabilities  with 
component  availabilities.  For  the  cases  in  which  the  com¬ 
ponents  do  not  operate  independently  and  each  component  has 
only  finitely  many  spares,  the  calculation  of  system  availa¬ 
bility  is  very  complex.  The  exact  expressions  depend  on  the 
nature  of  the  dependence.  We  explore  four  dependent  cases 
in  scenarios  2-5. 

Exact  expressions  for  system  availability  are  obtained  for 
the  two-component  series  module  described  by  scenarios  two 
and  three  when  there  are  infinitely  many  spares.  It  is  shown 
that  <_  A3<_  A2  and  (Y^J  _>  E[Y.j]  E[Y2J  where  Ai  is  the  steady 
state  availability  for  the  case  represented  by  scenario  i  and 
is  the  random  variable  representing  system  downtime.  For 
the  special  case  where  the  components  in  series  share  a  common 
pool  of  finitely  many  spares,  expressions  are  derived  for  the 
system  availability  in  all  five  scenarios.  The  results  obtained 
can  be  extended  to  series  systems  of  more  than  two  components 
with  proper  interpretation  of  the  scenarios. 

For  the  situation  in  which  each  component  has  its  own  unique 
pool  of  finitely  many  spares,  exact  computational  formulae  are 
derived  for  the  availability  of  series/parallel/mixed  system 
configuration  composed  of  independently  operating  components. 

For  scenarios  two  and  three  (symmetric  shut  down  and  asymmetric 
shut  down)  expressions  are  derived  for  the  Laplace  transforms 
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of  the  system  point  availability  for  any  given  vector  (n^,^) 
of  spares  allocation.  Numerical  examples  are  given  for  both 
scenarios  for  small  values  of  n^  and  nj*  For  larger  values 
of  n^  and  n ^  the  inversion  of  the  Laplace  transforms  is  diffi¬ 
cult.  However/  if  one  assumes  that  =  n 2  scenario  two, 
an  exact  expression  is  obtained  for  the  system  availability  for 
any  vector  (n^,n2) •  Other  approximations  for  system  availa¬ 
bility  are  provided  for  scenario  two  whenever  the  repair  rate 
is  much  larger  than  the  failure  rate.  A  numerical  example  of 
these  calculations  is  also  provided.  The  numerical  examples 
point  out  that  the  system  point  availability  is  much  higher 
when  components  share  a  common  pool  of  spares  than  when  the 
same  number  of  spares  is  available  but  no  sharing  is  allowed. 

For  scenarios  four  and  five  no  exact  expressions  are  determined 
for  system  availability.  However,  order  relationships  are 
obtained  which  allow  bounds  to  be  placed  on  system  availability. 

Because  some  of  the  expressions  for  system  availability 
under  the  different  scenarios  are  computationally  difficult 
it  might  be  useful  to  approximate  the  system  availability . 

Toward  that  objective,  order  relationships  among  the  scenario 
system  availabilities  are  obtained.  The  results  show  that 

AU»,m)  (t)  <  A^n'm)  (t)  <  A^n'm)(t) 


and 
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A^n,m)  (fc)  ^  A|n,m)  (t)  <  (t) 

r 

for  all  t  ^  0  and  n,m  e  IN  where  (n,m)  represents  the  vector 
of  spares  allocated  to  components  1  and  2  and  the  subscript 
refers  to  the  scenario.  The  same  order  relationships  hold 
when  components  share  a  common  pool  of  spares.  If  the  failure 
rates  are  small  compared  to  the  replacement  rates,  the  magnitude 
of  the  differences  is  small.  These  order  relationships  allow 
A|n,m)  A^n,m)  ^  to  uge(j  ag  boun(js  for  the  exact 

availability  for  the  other  scenario. 

The  final  product  of  this  research  is  an  algorithm  for 
allocating  resources  to  components  of  a  system  so  as  to  maxi¬ 
mize  system  availability  for  any  given  budget.  Kettelle  [221 
developed  an  algorithm  for  allocating  redundancy  to  the  com¬ 
ponents  of  a  series  system  where  all  components  in  series 
operate  independently.  The  algorithm  developed  in  this  paper 
utilized  Kettelle’s  idea  of  the  undominated  allocation  sequence 
and  his  way  of  tableau  computation  to  solve  the  more  general 
spares  allocation  problem  in  the  system  availability  optimi¬ 
zation  problem.  The  algorithm  groups  together  those  components 
which  depend  on  one  another  as  separate  independent  entities 
called  modules.  The  expressions  derived  in  this  thesis  are 
used  to  determine  the  availability  of  each  module  as  a  non¬ 
decreasing  function  of  the  resources  allocated  to  each  com¬ 
ponent  of  the  module.  The  algorithm  is  then  used  to  optimally 
allocate  resources  to  the  independent  modules.  The  resulting 


157 


algorithm  generalizes  Kettelle's  work  by  considering  system 
availability  vice  reliability  and  by  allowing  for  more  general 
series/parallel/mixed  system  configurations.  It  has  advan¬ 
tages  over  the  standard  dynamic  programming  algorithm  by 
eliminating  the  need  for  backtracking  and  by  solving  the 
allocation  problem  for  any  budget  <_  B.  By  careful  heuristic 
inspection  the  algorithm  can  be  made  very  efficient  for  manual 
calculations  because  large  blocks  of  cells  can  be  eliminated 
from  computation.  In  those  cases  where  two  stages  are  connected 
in  serires  the  remainder  of  an  entire  row  or  column  can  often 
be  eliminated  from  consideration.  Where  the  cumulative  stage 
returns  are  computed  from  parallel  connections  a  "north-west 
corner"  elimination  rule  can  sometimes  be  applied  to  eliminate 
entire  blocks  of  possibilities.  A  numerical  example  is  pro¬ 
vided  to  illustrate  the  allocation  algorithm.  The  example 
demonstrates  the  elimination  rules. 

The  exact  expressionsderived  in  this  study  are  sometimes 
computationally  intractable,  especially  for  automated  solution. 
In  actual  computations  some  of  the  terms  that  appear  in  the 
expressions  are  negligible.  Their  omission  simplifies  calcu¬ 
lation  dramatically.  Further  work  needs  to  be  done  in  deter¬ 
mining  decision  rules  which  will  indicate  when  an  expression 
can  be  ignored.  Also,  work  should  be  done  on  the  problem  of 
numerical  inversion  of  some  of  the  complicated  Laplace  trans¬ 
form  expressions .  Another  area  where  additional  research 
might  be  fruitful  is  the  development  of  a  computer  program  to 
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automate  the  allocation  algorithm.  A  computer  program  for 
implementing  the  algorithm  would  be  easy  to  develop  if  the 
program  is  to  ignore  the  computational  advantages  offered  by 
the  "elimination  rules."  A  program  which  systematically 
searches  for  an  opportunity  to  eliminate  blocks  of  possibili¬ 
ties  from  consideration  would  be  more  efficient;  albeit  more 
difficult  to  develop. 

This  research  points  out  that  system  availability  depends 
on  each  of  the  following  factors: 

(1)  Component  reliabilities 

(2)  Component  replacement  times 

(3)  System  configuration 

(4)  Operating  scenario  and  mission  profile 

(5)  Spare  parts  support  and  commonality  of  spare  parts 
These  factors  have  each  been  considered  explicitly  in  this 
thesis.  A  follow-on  study  should  be  made  to  determine  what 
data  are  needed  to  support  the  allocation  algorithm  developed 
in  this  thesis.  The  study  should  investigate  how  the  data 
should  be  collected.  A  data  base  management  system  should  be 
developed  to  allow  early  estimation  and  modification  of  and 
access  to  the  required  data.  An  important  reason  for  the 
simplistic  procedures  commonly  used  by  the  military  services 
in  constructing  allowance  lists  (see,  for  example,  the  U.S. 
Navy's  COSAL  or  FLSIP  models  [46])  is  because  the  data  requireed 
by  efforts  such  as  that  described  in  this  thesis  are  not  avail¬ 
able.  Any  model  or  algorithm  can  only  be  as  good  as  the  data 
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that  support  it.  Much  work  needs  to  be  done  in  this  area 
before  models  like  that  developed  here  can  be  implemented 
operationally. 

If  one  allows  for  maintenance  of  failed  components,  many 
additional  theoretically  challenging  problems  arise.  In 
actual  real-world  situations  it  may  be  that  a  component  is 
supported  not  only  by  spares,  but  also  by  a  maintenance  capa¬ 
bility,  where  malfunctioning  parts  are  repaired  and  held  for 
reuse.  The  addition  of  a  maintenance  capability  opens  up  a 
wealth  of  possibilities  by  allowing  for  different  maintenance 
policies,  number  of  servers,  service  discipline,  etc.  These 
are  proposed  as  fruitful  areas  of  future  study. 

Finally,  other  scenarios  describing  the  type  of  dependence 
that  may  exist  among  components  of  a  system  should  be  explored. 
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